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STATEMENT  OF  PROBLEM  STUDIED 


The  proposed  project  was  aimed  at  basic  understanding  of  the  nonlinear  force-displacement  and 
force-velocity  behavior  of  a  new  magneto-rheological  (MR)  fluid  damper.  The  major  research 
issues  and  the  objectives  of  this  project  were:  a)  to  develop  a  theoretical  study  to  predict 
nonlinear  behavior  of  new  MRF  dampers;  and  b)  to  conduct  a  comprehensive  experimental  study 
on  the  proposed  MRF  damper  to  validate  all  theoretical  findings.  Particularly,  the  focus  of  the 
research  was  on  the  specifications  required  for  the  Army’s  High  Mobility,  Multi-purpose 
Wheeled  Vehicle’s  (HMMWV’s)  suspension  system. 


SUMMARY  OF  IMPORTANT  RESULTS 


Introduction: 


Semi-active  vibration  control  devices  are  increasingly  being  investigated  and  implemented.  New 
advances  in  high-performance  materials  have  led  to  the  technological  growth  in  the  area  of  semi¬ 
active  protective  devices.  The  proposed  MRF  damper  system  will  take  advantage  of  the  state-of- 
the-art  in  advancement  of  materials  to  study  and  produce  a  practical,  efficient,  and  inexpensive 
shock  absorber  for  vibration  control  of  mechanical  systems. 

A  major  factor  in  tactical  doctrine  for  battlefield  operation  is  cross-country  mobility  for  combat 
vehicles  [1].  Advances  in  technology  that  allows  a  significant  enhancement  in  cross-country 
mobility  can  affect  the  battlefield  engagement’s  location  and  timing.  Extensive  test  results  by 
Army  TARDEC  have  clearly  demonstrated  that  the  utilization  of  semi-active  suspension 
technology  in  military  vehicles  can  increase  the  cross-country  mobility  more  than  30%  [2].  MR 
fluid  controllable  shock  absorbers  dramatically  reduce  the  mechanical  complexity  found  in 
current  servo-hydraulic  controllable  shock  absorbers.  Reservoirs,  heat  exchangers,  pumps  and 
external  hydraulic  systems  are  not  required  when  using  MR  technology.  The  fast  response 
behavior  of  MR  dampers  allow  for  active  control  throughout  a  higher  range  of  frequencies, 
therefore,  increasing  active  suspension  performance.  By  the  utilization  of  MR  fluid  dampers^ 
increased  reliability,  maintainability  and  performance  will  result.  Development  of  this  new  MR 
technology  for  vehicle  suspension  systems  can  matriculate  to  commercial  vehicles  as  well. 

In  order  to  accurately  define  the  design  specifications  of  the  new  University  of  Nevada,  Reno 
(UNR)  MR  fluid  damper,  Dr.  Francis  Hoogterp  of  the  U.S.  Army  TARDEC  was  consulted  for 
technical  data  and  performance  characteristics  of  HMMWV  suspension  components.  For  proper 
assessment  of  shock  absorber  performance  requirements,  two  original  equipment  manufacture 
(OEM),  Monroe  HMMWV  shock  absorbers  were  characterized  at  the  Composite  and  Intelligent 
Materials  Laboratory  (CIML,  http://web.me.unr.edu/ciml)  for  a  range  of  frequencies  and 
displacements.  These  test  results  provide  a  baseline  performance  standard  to  be  met  by  the  new 
MR  fluid  damper  design. 
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The  new  UNR  MR  damper  design  was  developed  to  meet  and  exceed  the  performance  criteria  set 
forth  by  the  OEM  test  results.  A  unique  capacity  for  bypass  valving  is  included  in  the  proposed 
design.  A  critical  element  of  vehicle  shock  absorbers,  bypass  valving  allows  the  shock  absorber 
to  accommodate  high-force  impulse  loading  (without  failure)  typical  to  off-highway 
environments.  The  new  MR  damper  has  bypass  valving  incorporated  into  the  design.  The  MR 
damper  was  tested  for  the  same  range  of  frequencies  and  displacements  as  the  OEM  units.  In 
addition,  the  same  tests  were  run  for  a  range  of  electrical  current  inputs,  the  means  of  controlling 
damping.  Moreover,  a  nonlinear  fluid-mechanics  based  theoretical  model  was  developed  by 
employing  Bingham  plastic  and  Herschel-Bulkey  non-Newtonian  fluid  models.  Three- 
dimensional  electromagnetic  finite  element  analyses  were  performed  for  establishing  a  base  for 
the  MR  fluid  damper  design. 

Several  prototype  UNR  MR  dampers  have  been  developed  and  tested  which  explore  the  design 
parameters.  These  dampers  were  tested  for  a  range  of  frequencies,  displacements  and  input 
currents.  Test  results  show  that  the  new  MR  fluid  shock  absorbers  can  provide  three  to  four 
times  the  amount  of  damping  force  than  provided  by  the  OEM  units.  In  addition,  without 
electrical  current  applied,  the  MR  fluid  shock  acts  as  a  passive  shock  absorber,  providing  the 
same  performance  as  the  OEM  unit. 

Because  the  MR  fluid  damper  uses  a  built-in  electromagnet  for  controlling  fluid  viscosity  (and 
ultimately  damping),  three-dimensional  magnetic  finite  element  analysis  was  performed  on  the 
device  to  determine  material  types,  magnetic  wire  size,  magnetic  saturation  of  fluid,  power 
requirements  and  optimal  geometry  for  the  best  MR-effect  (the  method  of  changing  fluid 
viscosity).  This  effort  began  with  development  of  a  three-dimensional  solid  model,  which  can  be 
directly  ported  into  different  finite  element  software  packages.  This  new  design  of  damper  has 
also  addressed  the  issue  of  fluid  sealing  by  adopting  a  seal  pack  that  does  not  leak  fluid.  This 
prototype  MR  fluid  damper  has  undergone  over  more  than  1000  different  tests  without  any  fluid 
leakage.  Also,  experimentally  investigated  is  the  use  of  a  satellite  accumulator.  Because  of  the 
nature  of  the  fluid,  correct  operating  pressure  for  an  MR  vehicle  shock  absorber  had  to  be 
determined  experimentally. 

This  summary  is  organized  in  two  sections: 

1 .  Nonlinear  theoretical  modeling  of  the  UNR  MRF  damper. 

2.  Experimental  studies. 


1.  THEORETICAL  MODELING: 


A  magneto-rheological  fluid  consists  of  micron-sized  ferrous  particles  within  a  base  carrier  fluid, 
such  as  silicon  or  mineral  oil,  and  additives.  Under  the  influence  of  a  magnetic  field,  the  ferrous 
particles  polarize  and  attract  each  other,  resulting  in  the  formation  of  chains  or  columns  of 
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particles  parallel  to  the  orientation  of  the  magnetic  flux  lines.  This  results  in  a  quantifiable 
change  in  the  “apparent  viscosity”  of  the  fluid,  which  is  defined  as: 


_  r 
Mapp  ~  T 

Y 

Here,  /uApp  is  the  apparent  viscosity  of  fluid,  x  represents  fluid  shear  stress  at  the  wall,  and 
is  the  shear  strain  rate.  The  Bingham  plastic  model  is  defined  as  [3]: 


(1) 
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where  is  the  field-induced  shear  stress  and  ^  represents  plastic  viscosity  of  the  fluid.  From 
Equations  (1)  and  (2),  one  has: 
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In  the  absence  of  minor  head  loss  terms,  as  well  as  fluid  inertia,  and  seal  friction,  the  pressure 
drop  across  the  piston  can  be  expressed  as  [4,5]0): 


ad _  ^>2fJ.AppLCHVCH  (4) 

V 

where  LCH  is  the  length  of  the  channel,  VCH  is  the  velocity  of  the  MRF  in  the  channel,  A#  is  the 
effective  hydraulic  diameter  (Deff  =  [64/(fRe)]D„ ,  where  JRe  is  the  geometry-specific,  laminar 
friction  constant  for  non-circular  flows),  and  A P  is  the  pressure  drop.  Substitution  of  Equation 
(3)  into  Equation  (4)  produces: 

aj}  _  ^2^LchVc„  |  22KIS Lm 

V  ^ 

where  Kf  is  the  shear  yield  stress  of  the  MRF  and  LMR  is  the  length  of  the  MR  valve.  I  is  the 
input  electric  current  and  K  and  Szxq  the  constants  which  are  determined  by  a  three-dimensional 
electro-magnetic  finite  element  analysis. 


For  different  rebound  and  compression  behavior,  or  to  accommodate  high-force  impact  loads, 
extra  passages  open  as  spring-backed  valves  retract.  For  this  MRF  damper  design,  these 
passages  are  in  a  parallel  arrangement.  Equating  total  volumetric  flow  rate  across  the  piston 
provides  the  following  relation  for  the  fluid  velocity  [6]: 


V  = 

VCH 


X_P  apiston 


NAmr  +  J3MAh 


where  is  the  input  velocity  of  piston,  Apjston  is  the  effective  area  of  piston,  Am  is  the  cross- 
sectional  area  of  MRF  channel,  Anon-mr  is  the  cross-sectional  area  of  non-MR-effected  passages, 
N  is  the  number  of  MR  channels,  M  is  number  of  non-MR-effected  passages  in  parallel,  /?  is  the 
contribution  (in  terms  of  percent-opening)  of  the  non-MR-effected  passages  and  ^  is: 


4 


(7) 


<*  = 


f  D1'" 

tzA 
D  2 

\u*ff  A 


I  I  £ ~yJy_  |  jAffl 

Ao  4 


(rr  /  /  A 

MqLp 


Here  D5  is  the  diameter  of  return  channels,  and  Lp  is  the  length  of  the  piston.  Combining 
Equations  (5)-(7)  and  multiplying  by  the  area  of  piston  yields  the  following  expression  for  the 
damping  force,  FDAmper‘ 
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Previous  studies  had  employed  the  Bingham  plastic  model.  A  non-dimensional  equation  of  flow 
was  introduced  by  Phillips  [3]  and  an  approximated  form  of  the  Phillips’s  theory  was  developed 
by  Gavin,  et  al  [7].  Because  of  its  simplicity,  the  Bingham  plastic  model  has  been  widely  used 
for  field-controllable  fluids,  however,  it  can  only  be  used  in  the  limiting  case  of  post-yield  phase 
when  the  material  exhibits  Newtonian  behavior.  ER  and  MR  fluids  are  generally  observed  to 
have  a  strong  field-dependent  shear  modulus  and  a  yield  stress  that  resists  the  material’s  flow 
until  shear  stress  reaches  a  critical  value.  The  Bingham  plastic  model  is  often  used  to  describe 
this  phenomenon.  In  the  Bingham  plastic  model,  ER  and  MR  fluids  are  assumed  to  be 
Newtonian  fluids  in  post-yield  regime,  with  a  constant  plastic  viscosity  assumption.  However, 
for  cases  where  the  fluid  experiences  post-yield  shear  thinning  or  shear  thickening,  the 
assumption  of  constant  plastic  viscosity  may  not  be  valid. 


In  this  paper,  the  previous  study  is  extended  to  include  generalized  analyses  of  ER  and  MR  fluids 
flowing  through  pipes  and  parallel  plates  using  the  Herschel-Bulkley  model.  Non-dimensional 
equations  are  derived,  and  simplified  closed-form  expressions  are  presented  for  determining  the 
pressure  drop  as  a  function  of  material  properties,  geometry,  and  volumetric  flow  rate.  The 
simplified  closed-form  Herschel-Bulkley  model  can  be  reduced  to  the  Bingham  plastic  model  for 
cases  where  post-yield  shear  thinning  or  thickening  are  minimal.  The  theoretical  Herschel- 
Bulkley  formulation  for  flow  through  pipes  is  compared  with  the  experimental  results  using  a 
MR  fluid  damper  which  is  constructed  and  tested  at  UNR.  Excellent  agreement  is  obtained. 


The  Herschel-Bulkley  constitutive  equation  presented  in  Equation  (1)  describes  the  flow  of  ER 
and  MR  fluids  through  a  circular  pipe  with  a  constant  cross  sectional  area. 


=  ry+k  | 
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—  =  0 
dr 
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dr 


\Trz\^T 
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(9) 


where  x n  is  shear  stress;  is  shear  strain  rate,  and  k  and  n  are  fluid  index  parameters.  t  v  is 

the  fluid  yield  stress  and  is  a  function  of  the  external  field.  Let  us  consider  a  steady,  one- 
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dimensional  flow  of  an  incompressible  ER  or  MR  fluid  through  a  straight  cylindrical  pipe  of 
constant  cross  section.  A  cylindrical  coordinate  system,  with  the  notation  shown  in  Figure  1,  is 

assumed.  The  momentum  equation  of  laminar  flow  in  a  pipe  can  be  written  as: 


rdp=  djrr J 
dz  dr 


where  the  pressure  gradient,  -j-  is  constant  along  the  flow  direction.  Integrating  Equation  (10) 
yields  shear  stress  as  a  function  of  radial  position: 


r„,_ir±+ a. 

2  dz  r 


The  axial  flow  in  the  pipe  requires  that  c0  be  zero.  The  flow  has  a  non-yield  region,  which  is 
signified  by  the  “plug”  radius  shown  in  Figure  1 .  The  plug  radius,  Rp,  can  be  determined  by: 


2  T 

R  =  y 
p  dp/dz 


(11a) 


The  yield  flow  region,  therefore,  is  defined  when  Rp  <r<R.  Substitution  of  Equation  (9)  into 


Equation  (11)  yields: 


T  +k±  „.Lr± 

dr  2  dz 


Since  — —  <  0 ,  Equation  (12)  can  be  written  as: 
dr 


1  „  dp  /  dz  O" 

2  k  k  ) 


Solving  Equation  (13)  with  boundary  condition  of  u  =  0  at r  =  R,  one  obtains  the  velocity 
distnbution  in  the  yield  flow  area  (R<r<R),  as  follows: 


u  _  n  2k  (  dp/dz  ty 
n  +  \  dp/dz  l  2k  V  k 


Ty  n 


dp/dz  R  r  ^  n 
2k  k  , 


In  the  plug  area,  0<r<Rp,  the  velocity  is  constant.  By  letting  r  =  Rp,  the  plug  velocity  can  be 
expressed  as: 
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The  volumetric  flow  rate  can  be  derived  from  the  following  relation: 


Q-  27t£  urdr  +  Rp2mip 


Substitution  of  Equations  (14)  and  (15)  into  Equation  (16)  yields: 
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where  p'=  -dp  /  dz .  Now,  let  us  consider  the  following  dimensionless  parameters: 

X  =  ^JL  =  ^-  and,  y  =  f— 1  ~  T  — 

p'D  R  \  n  nD* )  p'D 


(17a) 


Where  D—2R  and  X  represents  the  dimensionless  plug  thickness  or  dimensionless  shear  stress. 
When  X=l,  the  plug  radius  is  the  same  as  the  pipe  radius  and  the  yield  stress  of  the  fluid  is  equal 
to  its  shear  stress  at  the  wall  which  implies  that  there  is  no  flow  in  the  pipe.  Y  is  the 
dimensionless  volumetric  flow  rate.  By  using  the  above  dimensionless  variables,  Equation  (17) 
can  be  simplified  to  [7]: 


Y  =  (\-X)(\-aX-bX2 -cX3)n 


Where  a  = 
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Expanding  the  right  hand  side  of  Equation  (10)  by  Taylor  Series,  one  has: 
7=1-  (na  +l)X  +\  an  -  bn  +  2  )x  2  + 

\  2  J 
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For  X  <  0.5 ,  let  us  assume  a  linear  relationship  between  X  and  Y  in  the  form  of  1-Y=AX.  By 
linear  regression  of  Equation  (19),  the  coefficient  A  can  be  determined  as: 
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(19a) 


A  =  (na  +  \)--(na-bn+ — — —a2) 
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Therefore,  a  linear  approximation  to  Equation  (18),  for  all  n,  when  Y<  0.5,  can  be  presented  as: 
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Equation  (20)  represents  the  Herschel-Bulkley  fluid  model.  If  n=l,  the  Herschel-Bulkley  fluid 
equation  is  simplified  to  the  Bingham  plastic  fluid  model,  as  follows: 


7  =  1 — X  where  X  <  0.5 
3 


Equation  (20)  can  further  be  simplified  to  the  power  law  fluid.  By  setting,  t  =  0  (i.e.,  X=0), 

Equation  (20)  collapses  to  Y=l,  which  is  the  dimensionless  form  for  the  power  law  fluid  model. 
The  dimensioned  form  can  be  written  as: 


dp  _Ak(  8 Q  Yfjw  +  lY 
dz  D  v^£)3  J  V  n  j 

For  a  Newtonian  fluid,  if  n  is  set  to  1  in  Equation  (21a),  one  has: 

dp  _  128 kQ 
dz  kDa 


(21a) 


(21b) 


which  is  the  Hagen-Poiseuille  equation. 

Returning  to  the  Herschel-Bulkley  pipe  flow  for  ER  and  MR  field-controllable  fluids,  the 
dimensioned  form  of  Equation  (20)  provides  a  closed-form  expression  for  the  pressure  gradient 
as  follows: 


;  i<0.5 

dz  D  \  n  jdy )  D  R 

where^=3^±i.±i3^,l)(l-n) 

2»  +  l  16  (2n  +  l)2(«  +  l) 


(22) 


Equation  (22)  describes  the  pressure  loss  of  a  non-Newtonian  fluid  with  yield  stress.  The 
pressure  loss  can  be  separated  in  two  parts:  one  is  induced  by  pure  viscous  flow  and  the  other  is 
caused  by  yield  stress  when  the  non-yield  field  is  not  dominant  in  the  entire  flow  area.  When 
X>0.5,  the  plug  radius  reaches  half  of  the  pipe  radius,  the  yield  stress  will  play  an  important  role 
in  determining  pressure  loss.  A  simplified  expression  for  Equation  (18)  is  obtained  by 
considering  the  following  nonlinear  relation  between  X  and  Y: 
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Y  =  A'X(  i-xy 


(23) 


Where  A  and  a  are  functions  of  n  and  their  numerical  values  are  presented  in  Table  1.  The 
dimensioned  version  of  the  pressure  loss  can  be  derived  from  Equation  (23) 


dp 

dz 


where  A0  =  (A1)  “ 


D\ 


l -A 


4ry 

f  3n  + 1  8Q  V 

V  n  7d)i) 

0.5 


(24) 


The  exact.  Equation  (18),  and  the  simplified,  Equations  (20)  and  (23),  solutions  for  the  Herschel- 
Bulkley  pipe  flow  relationships  are  presented  in  Figure  2.  As  shown,  the  simplified  model  is  an 
excellent  approximation  to  the  exact  solution.  The  relative  maximum  error  between  the  two 
models  is  less  than  3%.  Therefore,  the  simplified  Herschel-Bulkley  relationships  presented  in 
quations  (22)  and  (24)  are  expected  to  provide  an  accurate  explicit  expression  for  the  pressure 
loss  associated  with  flow  of  ER/MR  fluid  in  pipes  with  circular  constant  cross  sections. 


2.  EXPERIMENTAL  STUDIES: 


To  establish  a  foundation  for  the  fail-safe  design,  first,  the  characteristics  of  two  OEM  HMMWV 
dampers  were  studied.  The  OEM  shock  absorbers  were  characterized  at  2,  3  4  and  7Hz 
frequencies.  For  each  frequency  considered,  four  different  peak-to-peak  amplitudes  of 
displacement  (0.5,  1.0, 1.5  and  2.0cm)  were  examined.  This  allowed  the  design  of  a  passive  fail¬ 
safe  baseline  characteristic  of  the  MRF  damper. 

For  the  MRF  damper,  an  additional  parameter  is  examined.  Increased  damping  is  achieved 
through  applying  electrical  current  to  the  built-in  electromagnet  of  the  device.  With  no  current 
supplied  to  the  damper  (i.e.,  zero-field,  or  passive-off),  the  MRF  damper  operates  in  fail-safe 
mode,  emulating  OEM  shock  absorber  behavior.  As  current  supplied  to  the  device  increases 
higher  damping  forces  result.  Three  different  electrical  current  inputs  are  investigated. 

Experiments  performed  on  a  shaking  table  (Ling  Dynamic  Systems,  model  V830-334T).  Figure 
3  is  an  illustration  of  the  experimental  test  setup.  A  linear  motion  transducer  (LMT)  was  used  to 
measure  displacement  and  velocity.  Electronic  signals  generated  by  the  load  cell  and  LMT  are 
sent  to  an  external  signal  processing  unit  which  scales  the  signal  for  the  data  acquisition  board  a 
National  Instruments  AT-MIO-16X  unit.  One  end  of  the  damper  is  mounted  directly  to  the 
shaking  table,  while  the  other  end  is  attached  to  a  13,000N  (3,0001b)  Eaton  Lebow  (Model  3132) 
load  cell.  A  personal  computer  houses  the  data  acquisition  board  and  the  instrumentation 
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software  package,  which  produce  test  data.  The  software  used  was  the  National  Instruments 
Lab  VIEW  software  package.  The  shake  table  is  controlled  by  another  computer,  which  controls 
power  and  frequency  to  the  shake  table  unit.  Power,  in  conjunction  with  frequency,  controls  the 
amplitude  of  table  travel,  which  is  monitored  by  the  LMT.  Equipment  were  calibrated  with 
known  masses,  frequencies,  measured  peak-to-peak  displacements  and  metered  electrical  current 
values. 

Two  identical  OEM  shock  absorbers  were  tested  to  assure  the  accuracy  of  test  results.  Each  test 
was  conducted  repetitively  for  a  duration  time  of  no  less  than  ten  seconds.  Figure  4  shows  an 
illustration  of  the  UNR  MRF  damper  along  side  (and  in  the  same  scale)  as  the  OEM  unit.  This 
illustration  shows  the  size  comparison  of  both  dampers.  The  geometric  constraints  of  the  shock 
absorber  include  length  of  travel,  mounting  means  and  outer  diameter  (inner  diameter  of  coil 
suspension  spring).  The  MRF  damper  was  constructed  specifically  for  testing  purposes  and  does 
not  have  all  of  the  elements  of  the  commercial  unit.  For  example,  elastomer  bearings  are 
incorporated  into  the  semi-rigid  mount  of  the  OEM  damper.  This  is  to  accommodate  for  minor 
off-axis  loading,  typical  to  pivotal  loading  environments  of  vehicles.  The  MRF  damper  has  rigid 
steel  mounts,  as  its  function  is  to  mount  only  to  the  testing  apparatus.  This  also  eliminates  the 
effects  of  the  spring  and  damping  properties  of  the  elastomer. 

Figure  5  is  a  photograph  of  the  actual  MRF  damper  used  for  testing.  Five  different  variants  of 
this  design  were  developed  and  tested  to  thoroughly  explore  parameters  of  design.  Figure  6  is  a 
photograph  of  the  satellite  accumulator,  which  is  kept  separate  from  the  experimental  damper  to 
simplify  damper  assembly  and  disassembly.  The  role  of  the  accumulator  is  to  accommodate  for 
the  change  in  available  fluid  volume  as  the  piston  reciprocates  inside  the  cylinder.  Most  OEM 
shock  absorbers  have  the  accumulator  integrated  into  the  unit  itself.  The  MRF  damper  design 
presented  in  this  paper  also  has  this  capacity,  but  the  accumulator  is  kept  separate  for  the 
aforementioned  purpose. 

Figure  7  presents  experimental  results  for  the  force-displacement  of  a  HMMWV  OEM  shock 
absorber  at  2Hz,  1  .Ocm  peak-to-peak  amplitude  of  displacement.  Figure  8  illustrates  the  force- 
velocity  data  for  the  same  test.  Figures  9  and  10  show  results  of  the  same  experiment  for  the 
UNR  HMMWV  MRF  damper.  In  comparing  Figure  5  and  Figure  9,  it  is  apparent  that  peak-to- 
peak,  passive-off  damping  force  is  similar  for  the  MRF  damper  and  the  OEM  unit.  Figure  9  also 
shows  the  range  of  damping  that  increases  as  current  supplied  to  the  MRF  damper  increases.  The 
baseline  peak-to-peak  force  at  2Hz  for  both  devices  is  about  600N  (1351bf),  while  the  increase  in 
peak-to-peak  force  (via  current  supplied  to  the  MRF  damper)  is  almost  2,200N  (5001bf).  In 
comparing  the  force- velocity  data  of  both  units  (Figure  6  for  the  OEM  unit  and  Figure  10  for  the 
MRF  damper  at  zero-field),  the  offset  ‘step’  in  the  hysteresis  loops  at  zero  velocity  is  due  to  the 
pressure  resulting  from  the  accumulator.  When  current  is  applied  to  the  MRF  damper,  damping 
force  increases  the  height  of  the  offset  force  step  at  zero  velocity.  In  Figure  10  the  peak  velocity 
increases  as  damping  force  increases.  This  is  the  result  of  the  experimental  setup.  The  shake 
table  maintains  constant  frequency,  while  the  peak-to-peak  displacement  changes,  slightly,  which 
results  in  increased  peak  velocity. 
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Figure  11  illustrates  another  means  for  evaluating  damper  performance.  Quantifying  the  area 
within  any  given  force-displacement  hysteresis  loop  provides  an  amount  of  work  done  by  the 
damper  during  one  complete  stroke.  For  an  equivalent  linear  damper,  the  damping  coefficient 
can  be  determined  by: 


,  W 

eg  nX2o)  (25) 


where  Ceq  is  the  equivalent  damping  coefficient,  W  represents  the  work  done  by  the  damper 
during  one  cycle,  X  is  the  amplitude  of  displacement,  and  co  is  the  frequency  in  radians  per 
second.  As  frequency  increases  and  current  to  the  damper  increases,  the  energy  absorbed  by  the 
damper  increases.  For  the  same  test  data,  Figure  12  illustrates  the  values  of  peak-to-peak  force  as 
a  function  of  input  electric  current  for  the  MRF  damper. 


Figure  13  and  Figure  14  illustrate  the  force  envelope  of  the  MRF  damper  as  a  function  of 
frequency  for  two  different  test  displacements  (1.0  and  1.5cm  peak-to-peak,  respectively).  Also 
shown  for  comparison  on  both  graphs  are  the  results  of  the  OEM  shock  absorber.  The  benefit  of 
controllable  damping  presented  in  this  study  is  made  evident  by  these  two  graphs;  a  range  of 
damping  force  with  the  MRF  damper  versus  the  passive  damping  force  of  the  OEM  shock 
absorber.  As  energy  inputs  vary  to  the  vehicle,  the  MRF  damper  has  the  capability  to  vary  the 
amount  of  energy  absorption  accordingly. 

Figure  15  is  a  comparison  of  force-displacement  experimental  and  theoretical  results  for  the 
MRF  damper  at  2Hz,  1.0cm  peak-to-peak  displacement.  Experiments  were  performed  at 
passive-off  mode  (zero  input  electrical  current)  and  3.0Amp  input  electric  current.  Figure  16 
shows  the  force-velocity  results  for  the  same  tests.  Figures  17  and  18  are  theoretical  force- 
displacement  results  illustrating  the  effects  of  the  extra  by-pass  valves  mentioned  earlier.  It  is 
evident  that  the  MRF  damper  no  longer  has  a  symmetrical  appearance  when  the  valves  are 
opened. 


Conclusions: 


New  MR  fluid  dampers  were  developed,  built,  tested  and  proven  effective.  The 
accomplishments  made  by  this  research  have  widespread  implications.  Commercial  vehicle 
suspension  systems,  including  automotive,  motorcycle  and  bicycle  applications  will  benefit  from 
an  MR  fluid  damper  design  that  has  incorporated  bypass  valving  into  its  design.  Once 
established  in  the  Army  HMMWV,  matriculation  of  an  MR-based  suspension  system  into 
civilian  applications  will  undoubtedly  result. 
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Figures  4.  MRP  (4a)  and  OEM  (4b)  dampers  at  the  same  scale. 


Figure  5.  Photograph  of  the  UNR 
MRF  HMMWV  damper. 


Figure  6.  Photograph  of  the  satellite 
accumulator. 
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Figure  7.  Force-displacement  loop  for  OEM  shock  absorber  at  2Hz,  1 ,0cm  p-t-p. 
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Figure  8.  Force- velocity  loop  for  the  OEM  shock  absorber  at  2Hz,  1 .0cm  p-t-p  displacement. 
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Figure  1 1 .  Experimental  results  for  equivalent  damping  coefficient  (Ceq)  versus  current  for 

different  frequencies  for  the  MRF  damper  tested  at  1.0cm  peak-to-peak  displacement. 
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Figure  12.  Experimental  results  for  peak-to-peak  force  as  a  function  of  input  electric  current  for  the 
MRP  damper  at  different  frequencies  (1.0cm  peak-to-peak  displacement). 
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Figure  14.  Comparison  of  peak-to-peak  force  as  a  function  of  frequency  for  the  OEM  unit  and  the 
MR  fluid  damper  at  1 .5cm  peak-to-peak  displacement. 
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Figure  15.  Comparison  between  theoretical  and  experimental  force  versus  displacement  of  the 
MRF  damper  for  passive-off  and  3.0Amp  input  electric  current. 
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Figure  16.  Comparison  between  theoretical  and  experimental  force  versus  velocity  of  the 
MRF  damper  for  passive-off  and  3.0Amp  input  current. 
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Figure  17.  Theoretical  force-displacement  results  for  the  MRP  damper  with  return  valves  opened  0.0% 
on  the  compression  stroke  at  3  Hz,  1.0cm  p-t-p,  for  passive-off  and  3.0Amp  electric  current. 
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Figure  1 8.  Theoretical  force-displacement  results  for  the  MRF  damper  with  return  valves  opened  25%  on 
the  compression  stroke  at  3Hz,  1.0cm  p-t-p,  for  passive-off  and  3.0Amp  electric  current. 
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ABSTRACT 

This  paper  presents  the  study  of  a  field-controllable,  semi-active  magneto-rheological  fluid 
(MRF)  shock  absorber  for  high-payload,  off-highway  vehicles.  A  MRF  damper  is  developed 
that  is  tailored  for  ground  vehicles  which  undergo  a  wide  range  of  dynamic  loading.  The  MRF 
damper  also  has  the  capability  for  different  rebound  and  compression  characteristics.  The  new 
MRF  shock  absorber  emulates  the  original  equipment  manufacture  shock  absorber  behavior  in  its 
passive-off  mode.  Theoretical  and  experimental  studies  are  performed  to  examine  this  MRF 
damper.  The  Bingham  Plastic  theory  is  employed  to  model  the  nonlinear  behavior  of  the  MRF. 
A  fluid-mechanics-based  theoretical  model  along  with  a  three-dimensional  finite  element 
electromagnetic  analysis  is  utilized  to  predict  the  MRF  damper  performance.  The  theoretical  and 
experimental  results  are  demonstrated  to  be  in  good  agreements. 


Keywords:  Fail-safe,  Field-controllable,  Semi-active,  Magneto-rheological  fluid  damper,  Off- 
highway,  Ground  vehicle,  HMMWV. 
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1.  INTRODUCTION 


Advances  in  technology  that  allow  a  significant  enhancement  in  cross-country  mobility  of 
combat  vehicles,  such  as  high  mobility  multi-purpose  wheeled  vehicles  (HMMWVs),  can  affect 
the  battlefield  engagement’s  location  and  timing  (Laiminie,  1992).  Experiments  have 
demonstrated  that  the  utilization  of  semi-active  suspension  technology  in  military  vehicles  can 
substantially  increase  the  cross-country  mobility  (Hoogterp  et  al.,  1993).  The  present  study 
provides  an  innovative,  fail-safe,  semi-active,  magneto-rheological  fluid  (MRF)  damper  for  the 
suspension  system  of  HMMWVs. 

To  achieve  HMMWV  s  high  performance  in  demanding  cross-country  terrain,  the  value  of  a 
properly  designed  and  reliable  suspension  system  cannot  be  over  estimated.  Conventional 
suspension  system  design  has  reached  its  practical  limits  towards  significantly  increasing  cross¬ 
country  performance  (Hoogterp  et  al.,  1993).  A  semi-active  suspension  system,  which  consists 
of  a  conventional  spring  and  a  controllable  damper,  has  shown  to  be  a  practical  solution  for 
enhancing  the  performance  of  suspension  systems.  The  variable  damper  provides  a  controllable 
damping  force.  It  has  also  been  demonstrated  that  semi-active  suspensions  can  provide  the  same 
performance  as  active  suspensions  (which  includes  power  actuators)  without  high  power 
consumption  and  with  very  small  design  variations  from  conventional  passive  systems  (Ivers, 
and  Miller,  1991,  and  Barak,  1989). 

As  a  component  of  a  semi-active  suspension  system,  the  role  of  the  MRF  shock  absorber  is  to 
have  the  capacity  to  produce  variable  damping  reactive  to  the  input  motion.  The  ability  to  vary 
damping  for  a  wide  range  of  input  motion  is  the  main  feature  of  semi-active  suspension  systems. 
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Other  features  include  control  of  body  roll  during  turning,  control  of  vehicle  pitch  during  braking 
and  acceleration,  and  accommodation  for  a  variable  vehicle  mass  (i.e.,  empty  or  full  payloads). 
In  comparison  to  typical  passive  suspension  systems,  semi-active  shock  absorbers  allow  for 
softer  damping  when  needed,  in  addition  to  harder  damping  for  situations  that  demand  it.  The 
main  benefits  of  this  technology  for  ground  vehicles  are:  faster  operational  speeds  over  similar 
terrain  and  improved  vibration  isolation  for  sensitive  payloads. 

MRF  dampers  are  semi-active  devices  that  contain  magneto-rheological  fluids  (Carlson,  et  al., 
1995).  The  activation  of  MRF  damper  s  built-in  magnetic  field  causes  a  fast  change  in  the  shear 
yield  stress  and  the  viscosity  of  a  MR  fluid.  The  fluid  changes  state  from  liquid  to  semi-solid  in 
milliseconds.  The  result  is  an  infinitely  variable,  controllable  damper  capable  of  large  damping 
forces.  MRF  dampers  offer  an  attractive  solution  to  energy  absorption  in  mechanical  systems 
and  structures.  This  is  because  they  can  be  battery  operated,  require  minimal  power  for 
operation,  and  have  additional  benefits,  such  as,  the  absence  of  mechanical  valving  (for  flow 
control)  in  the  damper,  insensitivity  to  impurities  penetration,  fluid  stability,  and  long  operational 
life.  Most  importantly,  MRF  dampers  can  be  considered  as  “fail-safe”  devices. 

This  study  presents  a  fail-safe  MRF  damper.  In  the  absence  of  an  electric  input  current  to  the 
damper  or  in  the  case  of  a  malfunction  in  the  control  hardware,  the  new  MRF  damper  provides 
the  same  damping  force  that  is  produced  by  an  original  equipment  manufacture  (OEM)  passive 
viscous  damper.  A  new  generation  of  MRF  dampers  has  been  developed  at  the  University  of 
Nevada,  Reno  (UNR)  that  provide  a  fail-safe  and  robust  vibration  damping  solution 
(Gordaninejad  and  Kelso,  1999,  and  Gordaninejad  and  Breese,  2000).  Several  innovative 
prototypes  MRF  dampers  have  been  developed  and  studied  in  recent  years  for  different 
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applications,  such  as  mountain  bicycles,  motorcycles,  vehicles,  and  seismic  applications  (Breese 
and  Gordaninejad,  1999,  Breese  and  Gordaninejad,  2000,  Breese,  et  aL,  2000,  Dogruoz,  et  al. 
2000,  Ericksen,  and  Gordaninejad,  1999,  Ericksen,  and  Gordaninejad,  2000,  Gordaninejad  and 
Breese,  2000,  Gordaninejad  and  Ericksen,  2000,  Gordaninejad,  and  Kelso,  2000,  Gordaninejad, 
et  al  1998  (a),  and  1998  (b),  Liu,  Y.,  et  al.,  Kelso  and  Gordaninejad,  1999,  Wang  and 
Gordaninejad,  2000  (a),  and  2000  (b)).  The  main  advantage  of  UNR  MRF  dampers  is  that  they 
can  produce  a  damping  force  nearly  the  same  as  an  OEM  passive  viscous  damper  at  their 
passive-off  state  (zero  electric  input  current)  while  providing  a  sizable  dynamic  force  range  to 
allow  controllability.  This  is  schematically  demonstrated  in  Figure  5,  where  the  OEM  viscous 
passive  damper  has  a  linear  response  and  the  MRF  damper  at  passive-off  state  is  designed  so  that 
it  is  slightly  softer  than  the  OEM  damper.  At  maximum  power  input,  the  damper  can  provide 
sufficient  dynamic  force  range  for  the  application. 

The  fail-safe  design  of  MRF  dampers  requires  that  the  device  must  have  a  certain  viscous 
damping  capability.  The  viscous  damping  force  in  a  damper  is  a  function  of  material  properties 
as  well  as  geometric  dimensions  within  the  damper.  This  requires  a  combined  investigation  of 
the  damper  and  the  fluid. 

In  this  work,  theoretical  and  experimental  studies  of  a  new  MRF  damper  are  presented 
(Gordaninejad  and  Kelso,  1999).  The  MRF  damper  can  provide  a  controllable  range  of  damping 
force  that  can  utilize  the  on-board  power  supply  of  the  HMMWV.  A  new  method  for  achieving 
controllable  damping  via  manipulation  of  MR  passage  geometry  allows  for  fine-tuning  of 
passive  and  semi-active  damping  behavior.  The  UNR  damper  design  also  has  the  capacity  to 
effectively  utilize  existing  ground  vehicle  shock  absorber  technology  in  the  arenas  of  by-pass 
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valves  and  over-travel  protection.  A  theoretical  model  is  developed  that  has  the  capability  of 
predicting  the  behavior  of  the  MRF  damper.  Utilizing  a  fluid-mechanics  based  model  and  the 
Bingham  plastic  theory,  the  theoretical  model  accommodates  for  the  geometric  parameters  as 
well  as  material  and  input  motion  characteristics.  In  addition,  the  theoretical  model  includes  the 
effects  of  bypass-valves. 


2.  THEORETICAL  MODELING  OF  THE  MRF  DAMPER 

A  magneto-rheological  fluid  consists  of  micron-sized  ferrous  particles  within  a  base  carrier 
fluid,  such  as  silicon  or  mineral  oil,  and  additives.  Under  the  influence  of  a  magnetic  field,  the 
ferrous  particles  polarize  and  attract  each  other,  resulting  in  the  formation  of  chains  or  columns 
of  particles  parallel  to  the  orientation  of  the  magnetic  flux  lines.  This  results  in  a  quantifiable 
change  in  the  “apparent  viscosity”  of  the  fluid,  which  is  defined  as: 


MaPP  —  . 

r 


(i) 


Here,  jiapp  is  the  apparent  viscosity  of  fluid,  t  represents  fluid  shear  stress  at  the  wall,  and  y 
is  the  shear  strain  rate.  The  Bingham  plastic  model  is  defined  as  (Phillips,  1969): 

t  =  Ty  +  /u0y  (2) 


where  xy  is  the  field-induced  shear  stress  and  jj$  represents  plastic  viscosity  of  the  fluid.  From 
Equations  (1)  and  (2),  one  has: 


In  the  absence  of  minor  head  loss  terms,  as  well  as  fluid  inertia,  and  seal  friction,  the  pressure 
drop  across  the  piston  can  be  expressed  as  (Kelso,  1998,  Kelso  and  Gordaninejad,  1999, 
Gordaninejad  and  Kelso,  2000): 


AP  = 


^Papp^ch  P CH 


where  LCh  is  the  length  of  the  channel,  VCH  is  the  velocity  of  the  MRF  in  the  channel,  A#  is  the 
effective  hydraulic  diameter  {Deff  —  [64/(fRe)]Dn ,  where  fRe  is  the  geometry-specific,  laminar 
friction  constant  for  non-circular  flows),  and  AP  is  the  pressure  drop.  Substitution  of  Equation 


(3)  into  Equation  (4)  produces: 


V  D«r 


where  Kl5  is  the  shear  yield  stress  of  the  MRF  and  Lmr  is  the  length  of  the  MR  valve.  /  is  the 
input  electric  current  and  K  and  5  are  the  constants  which  are  determined  by  a  three-dimensional 
electro-magnetic  finite  element  analysis. 


For  different  rebound  and  compression  behavior,  or  to  accommodate  high-force  impact  loads, 
extra  passages  open  as  spring-backed  valves  retract.  For  this  MRF  damper  design,  these 
passages  are  in  a  parallel  arrangement.  Equating  total  volumetric  flow  rate  across  the  piston 
provides  the  following  relation  for  the  fluid  velocity  (Kelso,  1998): 


V  = 

Y  CH 
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where  ^  is  the  input  velocity  of  piston,  AP1ST0N  is  the  effective  area  of  piston,  Am  is  the  cross- 
sectional  area  of  MRF  channel,  A non-mr  is  the  cross-sectional  area  of  non-MR-effected  passages, 
N  is  the  number  of  MR  channels,  M  is  number  of  non-MR-effected  passages  in  parallel,  /?  is  the 
contribution  (in  terms  of  percent-opening)  of  the  non-MR-effected  passages  and  $  is: 
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Here  Dp  is  the  diameter  of  return  channels,  and  Lp  is  the  length  of  the  piston.  Combining 
Equations  (5)-(7)  and  multiplying  by  the  area  of  piston  yields  the  following  expression  for  the 
damping  force,  F. damper' 
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3.  EXPERIMENTAL  STUDY 

To  establish  a  foundation  for  the  fail-safe  design,  first,  the  characteristics  of  two  OEM  HMMWV 
dampers  were  studied.  The  OEM  shock  absorbers  were  characterized  at  2,  3,  4  and  7Hz 
frequencies.  For  each  frequency  considered,  four  different  peak-to-peak  amplitudes  of 
displacement  (0.5,  1.0,  1.5  and  2.0cm)  were  examined.  This  allowed  the  design  of  a  passive  fail¬ 
safe  baseline  characteristic  of  the  MRF  damper. 

For  the  MRF  damper,  an  additional  parameter  is  examined.  Increased  damping  is  achieved 
through  applying  electrical  current  to  the  built-in  electromagnet  of  the  device.  With  no  current 
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not  have  all  of  the  elements  of  the  commercial  unit.  For  example,  elastomer  bearings  are 
incorporated  into  the  semi-rigid  mount  of  the  OEM  damper.  This  is  to  accommodate  for  minor 
off-axis  loading,  typical  to  pivotal  loading  environments  of  vehicles.  The  MRF  damper  has  rigid 
steel  mounts,  as  its  function  is  to  mount  only  to  the  testing  apparatus.  This  also  eliminates  the 
effects  of  the  spring  and  damping  properties  of  the  elastomer. 

Figure  3  is  a  photograph  of  the  actual  MRF  damper  used  for  testing.  Five  different  variants  of 
this  design  were  developed  and  tested  to  thoroughly  explore  parameters  of  design.  Figure  4  is  a 
photograph  of  the  satellite  accumulator,  which  is  kept  separate  from  the  experimental  damper  to 
simplify  damper  assembly  and  disassembly.  The  role  of  the  accumulator  is  to  accommodate  for 
the  change  in  available  fluid  volume  as  the  piston  reciprocates  inside  the  cylinder.  Most  OEM 
shock  absorbers  have  the  accumulator  integrated  into  the  unit  itself.  The  MRF  damper  design 

presented  in  this  paper  also  has  this  capacity,  but  the  accumulator  is  kept  separate  for  the 
aforementioned  purpose. 

4.  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

Figure  5  presents  experimental  results  for  the  force-displacement  of  a  HMMWV  OEM  shock 
absorber  at  2Hz,  1.0cm  peak-to-peak  amplitude  of  displacement.  Figure  6  illustrates  the  force- 
velocity  data  for  the  same  test.  Figures  7  and  8  show  results  of  the  same  experiment  for  the  UNR 
HMMWV  MRF  damper.  In  comparing  Figure  5  and  Figure  7,  it  is  apparent  that  peak-to-peak. 
passive-off  damping  force  is  similar  for  the  MRF  damper  and  the  OEM  unit.  Figure  7  also 
shows  the  range  of  damping  that  increases  as  cuirent  supplied  to  the  MRF  damper  increases. 
The  baseline  peak-to-peak  force  at  2Hz  for  both  devices  is  about  600N  (1351bf),  while  the 
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increase  in  peak-to-peak  force  (via  current  supplied  to  the  MRF  clamper)  is  almost  2,200N 
(5001bf).  In  comparing  the  force-velocity  data  of  both  units  (Figure  6  for  the  OEM  unit  and 
Figure  8  for  the  MRF  damper  at  zero-field),  the  offset  ‘step’  in  the  hysteresis  loops  at  zero 
velocity  is  due  to  the  pressure  resulting  from  the  accumulator.  When  current  is  applied  to  the 
MRF  damper,  damping  force  increases  the  height  of  the  offset  force  step  at  zero  velocity.  In 
Figure  8  the  peak  velocity  increases  as  damping  force  increases.  This  is  the  result  of  the 
experimental  setup.  The  shake  table  maintains  constant  frequency,  while  the  peak-to-peak 
displacement  changes,  slightly,  which  results  in  increased  peak  velocity. 


Figure  9  illustrates  another  means  for  evaluating  damper  performance.  Quantifying  the  area 

within  any  given  force-displacement  hysteresis  loop  provides  an  amount  of  work  done  by  the 

damper  during  one  complete  stroke.  For  an  equivalent  linear  damper,  the  damping  coefficient 
can  be  determined  by: 


kK2co  (9) 

where  Ceq  is  the  equivalent  damping  coefficient,  W  represents  the  work  done  by  the  damper 
during  one  cycle,  X  is  the  amplitude  of  displacement,  and  <o  is  the  frequency  in  radians  per 
second.  As  frequency  increases  and  current  to  the  damper  increases,  the  energy  absorbed  by  the 
damper  increases.  For  the  same  test  data,  Figure  10  illustrates  the  values  of  peak-to-peak  force 
as  a  function  of  input  electric  current  for  the  MRF  damper. 

Figure  11  and  Figure  12  illustrate  the  force  envelope  of  the  MRF  damper  as  a  function  of 
frequency  for  two  different  test  displacements  (1.0  and  1.5cm  peak-to-peak,  respectively).  Also 
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shown  for  comparison  on  both  graphs  are  the  results  of  the  OEM  shock  absorber.  Hie  benefit  of 
controllable  damping  presented  in  this  study  is  made  evident  by  these  two  graphs;  a  range  of 
damping  force  with  the  MRF  damper  versus  the  passive  damping  force  of  the  OEM  shock 

absorber.  As  energy  inputs  vary  to  the  vehicle,  the  MRF  damper  has  the  capability  to  vary  the 
amount  of  energy  absorption  accordingly. 

Figure  13  is  a  comparison  of  force-displacement  experimental  and  theoretical  results  for  the 

MRF  damper  at  2Hz,  1.0cm  peak-to-peak  displacement.  Experiments  were  performed  at 

passive-ofr  mode  (zero  input  electrical  current)  and  3.0Amp  input  electric  current.  Figure  14 

shows  the  force-velocity  results  for  the  same  tests.  Figures  15  and  16  are  theoretical  force- 

displacement  results  illustrating  the  effects  of  the  extra  by-pass  valves  mentioned  earlier.  It  is 

evident  that  the  MRF  damper  no  longer  has  a  symmetrical  appearance  when  the  valves  are 
opened. 

5.  CONCLUSIONS 

A  new,  fail-safe,  controllable,  semi-active  MRF  fluid  damper  was  developed  and  studied  to 
determine  the  viability  of  implementation  of  this  technology  into  ground  vehicles.  It  has  been 
demonstrated  that  the  MRF  damper  has  a  broad  dynamic  force  range  and  can  emulate  OEM 
damper  behavior  at  its  passive-off  fail-safe  mode.  The  new  design  also  allows  for  integration  of 
cmrent  hydraulic  automotive  shock  absorber  technology,  such  as  by-pass  valves  and  over-travel 

protection.  In  conjunction  with  a  semi-active  control  system,  this  technology  has  shown  promise 
for  implementation  into  ground  vehicles. 
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Figure  1.  Schematic  of  the  experimental  setup, 


Figures  2.  MRF  (2a)  and  OEM  (2b)  dampers  at  the  same  scale. 


Figure  3 .  Photograph  of  the  UNR 
MRF  HMMWV  damper, 


Figure  4.  Photograph  of  the  satellite 
accumulator. 
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Figure  7 .  Force-displacement  loops  for  the  MRF  damper  at  2Hz,  1  .Ocm  p-t-p 
displacement  for  different  electrical  current  inputs. 
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Figure  8.  Force-velocity  loops  for  MRF  damper  at  2Hz,  1.0cm  p-t-p  displacement  for 
different  electrical  current  inputs. 
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Experimental  results  for  equivalent  damping  coefficient  (C^)  versus  current  for  different 
frequencies  for  the  MRF  damper  tested  at  1.0cm  peak-to-peak  displacement. 
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Figure  1 0.  Experimental  results  for  peak-to-peak  force  as  a  function  of  input  electric  current  for  the 
MRF  damper  at  different  frequencies  (1 .0cm  peak-to-peak  displacement). 
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Figure  1 1 .  Comparison  of  peak-to-peak  force  as  a  function  of  frequency  for  the  OEM  unit  and  the 
MR  fluid  damper  at  1  .Ocm  peak-to-peak  displacement. 


'  X--  1  —  -  j  -  ;  7~~~  '• 


'  ...  r.  ..  ‘ 


i*  *  ~ 


Dynamic 

Force 

Range 


— O  .MRN-1-2,  O.OAmp 
-  MRN-1-2, 3.0Amp 


Test  Frequency  (Hz) 


Figure  12.  Comparison  of  peak-to-peak  force  as  a  function  of  frequency  for  the  OEM  unit  and  the 
MR  fluid  damper  at  1 ,5cm  peak-to-peak  displacement. 
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Figure  13.  Comparison  between  theoretical  and  experimental  force  versus  displacement  of  the 
MRF  damper  for  passive-off  and  3.0Amp  input  electric  current. 


1 - - - - - .3000  _ _ _ _ _ _ 

Velocity  (cm/s) 

Figure  14.  Comparison  between  theoretical  and  experimental  force  versus  velocity  of  the 
MRF  damper  for  passive-off  and  3.0Amp  input  current. 
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Figure  15.  Theoretical  force-displacement  results  for  the  MRF  damper  with  return  valves  opened  0.0% 
on  the  compression  stroke  at  3  Hz,  1.0cm  p-t-p,  for  passive-off  and  3.0 Amp  electric  current. 
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Figure  16.  Theoretical  force-displacement  results  for  the  MRF  damper  with  return  valves  opened  25%  on 
the  compression  stroke  at  3Hz,  1.0cm  p-t-p,  for  passive-off  and  3.0Amp  electric  current. 
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ABSTRACT:  The  Bingham  plastic  constitutive  model  has  been  widely  used  to  predict  the 
post-yield  behavior  of  electro-  and  magneto-rheological  fluids  (ER  and  MR  fluids).  However,  if 
these  fluids  experience  shear  thinning  or  shear  thickening,  the  Bingham  plastic  model  may  not  be  an 
accurate  predictor  of  behavior,  since  the  post-yield  plastic  viscosity  is  assumed  to  be  constant.  In  a  re¬ 
cent  study,  it  was  theoretically  and  experimentally  demonstrated  that  the  Herschel-Bulkley  fluid 
model  can  be  successfully  employed  when  evaluating  non-Newtonian  post-yield  behavior  of  ER  and 
MR  fluids.  In  this  paper,  we  extend  our  previous  work  and  adopt  the  Herschel-Bulkley  model  to  in¬ 
clude  a  detailed  analysis  of  ER  and  MR  fluid  dynamics  through  pipes  and  parallel  plates.  Simplified 
explicit  expressions  for  the  exact  formulation  are  also  developed.  It  is  shown  that  the  proposed  sim¬ 
plified  model  of  the  Herschel-Bulkley  steady  flow  equations  for  pipes  and  parallel  plates  can  be  used 
as  an  accurate  design  tool  while  providing  a  convenient  and  generalized  mathematical  form  for  mod¬ 
eling  ER  and  MR  fluids.  Theoretical  and  experimental  analyses  are  presented  for  a  MR  fluid  damper, 
which  is  designed,  developed,  and  tested  at  the  University  of  Nevada.  Reno  (UNR). 


INTRODUCTION 

Electro-  and  magneto-rheological  (ER  and  MR)  fluids 
demonstrate  non-linear  material  behavior  when  sub¬ 
jected  to  an  external  electric  or  magnetic  field,  respectively. 
The  rheological  behavior  of  these  materials  can  be  separated 
into  two  distinct  pre-  and  post-yield  regimes.  In  most  de¬ 
vices,  such  as,  dampers,  valves,  and  clutches,  the  post-yield 
behavior  of  ER  and  MR  fluids  is  the  dominant  mode  of  oper¬ 
ation.  When  the  fluid  undergoes  a  shear  thinning  or  thicken¬ 
ing,  its  post-yield  behavior  is  non-linear.  This  effect  requires 
a  fluid  model  to  account  for  the  non-Newtonian  behavior  of 
ER  and  MR  fluids  in  the  post-yield  regime. 

Previous  studies  have  employed  the  Bingham  plastic 
model.  A  non-dimensional  equation  of  flow  was  introduced 
by  Phillips  (1969)  and  an  approximated  form  of  the  Phillips’s 
theory  was  developed  by  Gavin,  et  al.  (1996).  Because  of  its 
simplicity,  the  Bingham  plastic  model  has  been  widely  used 
for  field-controllable  fluids,  however,  it  can  only  be  used  in 
the  limiting  case  of  post-yield  phase  when  the  material  ex¬ 
hibits  Newtonian  behavior.  ER  and  MR  fluids  are  generally 
observed  to  have  a  strong  field-dependent  shear  modulus  and 
a  yield  stress  that  resists  the  material’s  flow  until  shear  stress 
reaches  a  critical  value.  The  Bingham  plastic  model  is  often 
used  to  describe  this  phenomenon  (Stanway  et  al.,  1996).  In 
the  Bingham  plastic  model,  ER  and  MR  fluids  are  assumed  to 
be  Newtonian  fluids  in  post-yield  regime,  with  a  constant 
plastic  viscosity  assumption.  However,  for  cases  where  the 
fluid  experiences  post-yield  shear  thinning  or  shear  thicken- 
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ing,  the  assumption  of  constant  plastic  viscosity  may  not  be 
valid. 

Goodwin  et  al.  (1997)  reported  that  the  ER  fluids  prepared 
with  particles  of  low  conductivity  showed  pseudo-plastic  be¬ 
havior  with  the  power  law  index  that  was  a  function  of  con¬ 
ductivity.  Mokeev  et  al.  (1992)  employed  the  continuous 
inhomogeneous  dielectric  medium  approximation  to  model 
Couette  shear  in  high-concentration  ER  fluid  under  electric 
field.  They  obtained  an  exponential  dependence  of  the  effec¬ 
tive  viscosity  on  the  electric  intensity  and  shear  rate,  which  in 
the  limiting  case  is  reduced  into  the  Bingham  plastic  model. 
Shulman  and  Korobko  (1978)  considered  refining  the 
Bingham  plastic  model  into  a  four-parametric  Casson-type 
model  with  coefficients  depending  on  electric  intensity. 
Stanway  et  al.  (1989)  proposed  a  model  to  describe  control¬ 
lable  fluid  damper  behavior.  An  experimental  technique  of 
nonlinear  sequential  filtering  was  used  to  estimate  the  pa¬ 
rameters  associated  with  an  nth-power  velocity  model  for  the 
damping  mechanism. 

In  a  recent  work  by  Marksmeier  et  al.  (1998),  a  new  ER 
grease  material  was  evaluated  at  UNR.  A  prototype  ER 
grease  damper  was  developed  and  tested.  The  shear  thinning 
effect  of  ER  grease  was  more  pronounced  than  that  of  regular 
ER  fluid.  In  order  to  accurately  model  the  damper  behavior, 
the  Herschel-Bulkley  model  was  used.  This  model  allows  for 
a  non-linear  post-yield  behavior.  Subsequently,  by  using  a 
lumped  parameter  model  for  the  damper,  the  force-displace¬ 
ment  and  force-velocity  relations  were  established.  The 
comparisons  between  experimental  results  and  theoretical 
analyses  utilizing  both  the  Herschel-Bulkley  and  the 
Bingham  plastic  models  confirmed  that  the 
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Herschel-Bulkley  model  predictions  are  more  accurate  than 
the  Bingham  plastic  model.  It  was  also  concluded  that  the 
Herschel-Bulkley  fluid  model  is  a  generalized  model  for 
viscoplastic  flow  with  yield  stress  and  can  be  a  valuable  tool 
in  the  analysis  and  design  of  controllable  fluid  dampers. 

In  this  paper,  the  previous  study  is  extended  to  include 
generalized  analyses  of  ER  and  MR  fluids  flowing  through 
pipes  and  parallel  plates  using  the  Herschel-Bulkley  model. 
Non-dimensional  equations  are  derived,  and  simplified 
closed-form  expressions  are  presented  for  determining  the 
pressure  drop  as  a  function  of  material  properties,  geometry, 
and  volumetric  flow  rate.  The  simplified  closed-form 
Herschel-Bulkley  model  can  be  reduced  to  the  Bingham 
plastic  model  for  cases  where  post-yield  shear  thinning  or 
thickening  are  minimal.  The  theoretical  Herschel-Bulkley 
formulation  for  flow  through  pipes  is  compared  with  the  ex¬ 
perimental  results  using  a  MR  fluid  damper  which  is  con¬ 
structed  and  tested  at  UNR.  Excellent  agreement  is  obtained. 

MODELING  OF  ER  AND  MR  FLUID  FLOW  USING 
THE  HERSCHEL-BULKLEY  MODEL 

In  this  section,  formulations  for  pipes  and  parallel  plates 
are  presented.  In  addition,  a  simplified  model  is  developed 
for  each  case.  It  is  demonstrated  that  the  Herschel-Bulkley 
formulation  can  collapse  to  the  Bingham  plastic  model  for 
fluid  index  number  of  n  =  1 . 


Figure  1.  ER  or  MR  fluid  flow  profile  in  a  pipe  with  a  circular  constant 
cross  section. 


1  dp  Cn 

^=~2rJz+  7  & 

The  axial  flow  in  the  pipe  requires  that  c0  be  zero.  The  flow 
has  a  non-yield  region,  which  is  signified  by  the  “plug”  ra¬ 
dius  shown  in  Figure  1 .  The  plug  radius,  Rp ,  can  be  deter¬ 
mined  by: 


Flow  Through  a  Circular  Pipe  with  a 
Constant  Cross  Section 


Rp  dp/dz 


The  Herschel-Bulkley  constitutive  equation  presented  in 
Equation  (1)  describes  the  flow  of  ER  and  MR  fluids  through 
a  circular  pipe  with  a  constant  cross  sectional  area. 


The  yield  flow  region,  therefore,  is  defined  when  Rp<r< 
R.  Substitution  of  Equation  (1)  into  Equation  (3)  yields: 


,  uu 

Xr.=X,.+k —  x  >T 
*  •’  dr 


-=o 

dr 


No¬ 


where  Tr;  is  shear  stress;  du/dr  is  shear  strain  rate,  and  k  and  n 
are  fluid  index  parameters.  xv  is  the  fluid  yield  stress  and  is  a 
function  of  the  external  field.  Let  us  consider  a  steady, 
one-dimensional  flow  of  an  incompressible  ER  or  MR  fluid 
through  a  straight  cylindrical  pipe  of  constant  cross  section. 
A  cylindrical  coordinate  system,  with  the  notation  shown  in 
Figure  1,  is  assumed.  The  momentum  equation  of  laminar 
flow  in  a  pipe  can  be  written  as: 


rdp=  d(rt  ,,) 
dz  dr 


where  the  pressure  gradient,  dp/dz  is  constant  along  the  flow 
direction.  Integrating  Equation  (2)  yields  shear  stress  as  a 
function  of  radial  position: 


2  r  dz 


Since  du/dr  <  0,  Equation  (4)  can  be  written  as: 


2  dp/dz  Ty  » 
2 r  k  k  , 


Solving  Equation  (5)  with  boundary  condition  of  u  =  0 
at  r=R,  one  obtains  the  velocity  distribution  in  the  yield  flow 
area  (Rp  <r<R),  as  follows: 


n  2k 
n  + 1  dp/dz 


dp/dz  Tv  |  " 
2k  r~~k  , 


dp/dz  R  Tv  I  n 
2k  k  , 


In  the  plug  area,  0  <r<Rp,  the  velocity  is  constant.  By  let¬ 
ting  r  =  Rp,  the  plug  velocity  can  be  expressed  as: 
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n  2k 
n  + 1  dp!  dz 


dp/dz  h 

~*'~7 


dp/dz  _  % 

- A - 

2k  k 


For  X  <  0.5,  let  us  assume  a  linear  relationship  between  X 
and  Fin  the  form  of  1  -  Y=AX.  By  linear  regression  of  Equa¬ 
tion  (1 1),  the  coefficient  A  can  be  determined  as: 


A  =  (na  + 1)-—  na-bn+ — — —a2  (11a) 

81  2  ) 


The  volumetric  flow  rate  can  be  derived  from  the  follow¬ 
ing  relation: 


Q  =  2njR  urdr  +  RpKUp 


Substitution  of  Equations  (6)  and  (7)  into  Equation  (8) 
yields: 


(Pl.x)  2xf^-T 

O'R  )  n  3  - L  + - Li - 

— — Tv  n  R3  3/1  +  1  2/1  +  1 


Therefore,  a  linear  approximation  to  Equation  (10),  for  all 
n,  when  X  <  0.5,  can  be  presented  as: 


Y  (na  +  \)-~  na-bn+ — — — fl2  X  (12) 
81  2  ) 


Equation  (12)  represents  the  Herschel-Bulkley  fluid 
model.  If  n  =  1,  the  Herschel-Bulkley  fluid  equation  is  sim¬ 
plified  to  the  Bingham  plastic  fluid  model,  as  follows: 


Y  =  1  —  X  where  X<0.5 
3 


Equation  (12)  can  further  be  simplified  to  the  power  law 
fluid.  By  setting,  Tv  =  0  (i.e.,  X  =  0),  Equation  (12)  collapses 
to  Y  =  1,  which  is  the  dimensionless  form  for  the  power  law 
fluid  model.  The  dimensioned  form  can  be  written  as: 


where  p '  =  - dp/dz .  Now,  let  us  consider  the  following 
dimensionless  parameters: 


dp  _Ak(  8  Q 
dz  D\kD3 


«_  (9a) 

pD  R  {  n  kD3)  p'D 


For  a  Newtonian  fluid,  if  n  is  set  to  1  in  Equation  ( 1 3a),  one 
has: 


where  D-2R  and X represents  the  dimensionless  plug  thick¬ 
ness  or  dimensionless  shear  stress.  When  X=  1,  the  plug  ra¬ 
dius  is  the  same  as  the  pipe  radius  and  the  yield  stress  of  the 
fluid  is  equal  to  its  shear  stress  at  the  wall  which  implies  that 
there  is  no  flow  in  the  pipe.  Y is  the  dimensionless  volumetric 
flow  rate.  By  using  the  above  dimensionless  variables,  Equa¬ 
tion  (9)  can  be  simplified  to  (Chilton  and  Stainsby,  1998): 


Y  =  ( 1  -  X)(  1  -  aX  -  bX2  -  cX3 )" 


where 


f-JJm  (13b) 

dz  7i  D4 

which  is  the  Hagen-Poiseuille  equation. 

Returning  to  the  Herschel-Bulkley  pipe  flow  for  ER  and 
MR  field-controllable  fluids,  the  dimensioned  form  of  Equa¬ 
tion  (12)  provides  a  closed-form  expression  for  the  pressure 
gradient  as  follows: 


dp  ^  4tv  (ln  +  \  8Q  Y1 


n  nD 


— ;  —  <  0.5  (14) 


-  1  b  ~ 

a~2n  +  \'  ~  (n  +  l)(2/i  +  l)’ 


(n  +  l)(2n  + 1) 


Expanding  the  right  hand  side  of  Equation  (10)  by  Taylor 
Series,  one  has: 


Y  =  l-(na  +  l)X+fan-bn+— — —a2  X2 

V  2  y 

+  c(  (~nc+ bn +n(n- 1)  -  —  - 1  j-  -  -  a3 )  X3  +  -  (11) 
V  6  ) 


where 

3n  +  l  3  (3/7  + 1)(1  —  n) 

2/J  + 1  16  (2/7  +  l)2(rt  +  l) 

Equation  (14)  describes  the  pressure  loss  of  a  non-Newto¬ 
nian  fluid  with  yield  stress.  The  pressure  loss  can  be  sepa¬ 
rated  in  two  parts:  one  is  induced  by  pure  viscous  flow  and 
the  other  is  caused  by  yield  stress  when  the  non-yield  field  is 
not  dominant  in  the  entire  flow  area.  When  X  >  0.5,  the  plug 
radius  reaches  half  of  the  pipe  radius,  the  yield  stress  will 
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play  an  important  role  in  determining  pressure  loss.  A  sim¬ 
plified  expression  for  Equation  ( 1 0)  is  obtained  by  consider¬ 
ing  the  following  nonlinear  relation  between  X  and  Y: 

Y  =  A'X(l-X)a  (15) 


where  vv  is  the  width  of  parallel  plates  and  zp  is  the 
dimensionless  plug  thickness  defined  as: 


=  2x1  =  ^ 

^  p'h  h 


(17a) 


where  A'  and  a  are  functions  of  n  and  their  numerical  values 
are  presented  in  Table  1 .  The  dimensioned  version  of  the 
pressure  loss  can  be  derived  from  Equation  (15) 


Here,  h  is  the  distance  between  the  two  parallel  plates,  and  hp 
is  the  plug  thickness.  Let  us  define  the  following 
dimensionless  variable: 


dp 

dz 


4tv 


D 


1-4 


(3b +  1 

8  Q  V 

V 

l  » 

7tD3J 

for 


R 


>0.5 


(16) 


where 


V  = 


2/7  +  1 
In 


Y  2k_ 
\wh2  J  p'h 


(17b) 


Using  Equation  (17b),  Equation  (17)  can  be  rewritten  as: 


V  =  (l-zJ 


M+l 


n  +  l 


+  1 


(18) 


Ao=(A')  a 

The  exact,  Equation  (10),  and  the  simplified,  Equations 
(12)  and  (15),  solutions  for  the  Herschel-Bulkley  pipe  flow 
relationships  are  presented  in  Figure  2.  As  shown,  the  simpli¬ 
fied  model  is  an  excellent  approximation  to  the  exact  solu¬ 
tion.  The  relative  maximum  error  between  the  two  models  is 
less  than  3%.  Therefore,  the  simplified  Herschel-Bulkley  re¬ 
lationships  presented  in  Equations  (14)  and  ( 1 6)  are  expected 
to  provide  an  accurate  explicit  expression  for  the  pressure 
loss  associated  with  the  flow  of  ER/MR  fluid  in  pipes  with 
circular  constant  cross  sections. 

Flow  through  Parallel  Plates 

A  similar  approach  used  for  pipes  with  circular  cross  sec¬ 
tions  can  be  employed  to  develop  relations  for 
Herschel-Bulkley  flow  between  two  parallel  plates.  Let  us 
assume  that  the  parallel  plates  are  fixed  and  a  constant  pres¬ 
sure  gradient  p  provides  the  fluid  flow  between  the  plates. 
The  Herschel-Bulkley  solution  for  steady,  one-dimensional 
flow  with  no-slip  boundary  condition  yields  the  following 
volumetric  flow  rate: 

(1  —  ~„)2  t  2;,(1  -Zp)  |  4 
3  n  + 1  2n  +  l  n  + 1 
n  n  n  . 

(17) 


Expanding  Equation  (18),  one  has: 

V  =  1 +clz„+c2z2p+c3zl  -  (19) 

where 

2n  +  l  >i(l  +  2n)(l-n)  _  «(n2  +1) 

C|  _-7+T;  C2~  2(n  + 1)2  ’  C'  2(n  +  l) 

For  zp  <  0.5,  using  the  same  method  as  described  in  the  pre¬ 
vious  section,  the  linear  approximation  of  Equation  (19)  can 
be  expressed  as: 

v  =  1+(C|+iC2+ir3)v’  Zp~0'5  (20) 

For  zp  >  0.5,  a  similar  expression  to  the  one  given  in  Equa¬ 
tion  ( 1 5)  can  be  obtained,  as  follows: 

V  =  B'zp{\-zpf ;  zp  > 0.5  (21) 

where  B'  and  P  are  functions  of  power  index,  n,  that  are 
obtained  by  curve  fitting  to  analytical  solutions  of  Equation 
(18),  and  are  presented  in  Table  2.  Again,  similar  to  the 
pipe  flow  case.  Figure  3  demonstrates  the  simplified 
Herschel-Bulkley  model  for  laminar  flow  between  the  paral¬ 
lel  plates.  Equations  (20)  and  (21)  show  an  excellent  ap¬ 
proximation  to  the  exact  formulation  presented  in  Equa¬ 
tion  ( 1 8).  The  dimensioned  forms  of  Equations  (20)  and  (2 1 ) 
are: 


4 Q_ 
wh 2 


2k 

p'h 


n+ 1 

=  d -zp)n 


Table  1.  Values  of  A',  a  and  Ao  for  different  n. 


n 

a 

A 

Ao  =  (A')"1/a 

1 

2.260 

3.289 

0.590 

0.8 

2.105 

3.031 

0.591 

0.4 

1.784 

2.700 

0.573 

Table  2.  Values  of  B',  P  and  B0  for  different  n. 


n 

p 

B' 

B0  =  (B')-1,e 

1 

2.629 

3.952 

0.593 

0.8 

2.475 

3.810 

0.583 

0.4 

1.816 

2.727 

0.575 
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lure  2.  Comparison  between  the  exact  and  simplified  versions  of  Herschel-Bulkley  pipe  flow  formulation  fortwo  different  flow  indices,  n  - 
nd  n  =0.4. 


Figure  3.  Comparison  between  the  exact  and  simplified  versions  of  Herschel-Bulkley  parallel  plates  flow  formulation  fortwo  different  fluid  in- 
dices,  n  =  1  and  n  =0.4. 
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_ Herschel-Bukley  model  (k  =  430,  n  -  0.4) 


Figure  6.  Shear  stress  versus  shear  strain  rate  at  various  magnetic 
field  strengths. 


of  low  magnetic  field  strength,  MR  fluids  exhibit  Bingham 
plastic  behavior  with  a  constant  plastic  viscosity  being  equal 
to  the  zero-field  viscosity  of  Newtonian  fluids.  However,  at 
higher  magnetic  field  strengths,  the  MR  fluid  exhibits 
pseudo-plastic  behavior  with  a  field-dependent  yield  stress. 
Considering  the  shear  thinning  effect,  the  Herschel-Bulkley 
model  is  adapted  (where  k = 430  and  n = 0.4)  to  represent  the 
measured  shear  stress  versus  shear  strain  rate  data  (see  Fig¬ 
ure  6)  The  yield  stress  Xv  is  a  function  of  the  applied  mag¬ 
netic  field  strength.  For  this  UNR  MR  fluid  damper  the  rela¬ 
tionship  between  the  yield  stress  xv  and  the  input  current  I  is: 

tv  =  2.7xl04/15  (26) 

where  the  units  for  xv  and  /  are  Pa  and  amp,  respectively. 
Equation  (26)  is  conformed  with  both  the  experimental  and 
theoretical  results  that  the  yield  stress  x,  is  a  power-law  func¬ 
tion  of  the  flux  density,  xv «  B1 5  (Ginder,  1998). 

The  overall  pressure  drop  A P  across  the  piston  is  mainly 
contributed  by  two  parts,  one  is  the  pressure  drop  A Pris  which 
is  due  to  the  fluid  viscosity,  and  the  other  is  the  pressure  drop 
A p  aCross  the  MR  valve  normal  to  the  axis  of  the  shaft.  The 
totalr pressure  drop  across  the  piston  head  can  be  written  as 
following: 


AP  =  APvis+2APmr  (27) 


For  a  Newtonian  Poiseuille  flow  in  a  pipe,  the  pressure 
drop  APy-is  is  expressed  as: 


A Pr„  = 


128  pQA 


(28) 


The  MR  valve  pressure  drop  A Pmr  can  be  calculated  by 
Equations  (14)  or  (16),  which  depends  on  the  value  of 
nondimensional  parameter  X.  Here,  a  new  dimensionless  pa¬ 
rameter  is  introduced: 


T  =■ 


*v 


(3n  +  l)  8 Q 


V 


(29a) 


n  nD 

Since  t  =  XI Y,  then,  from  Equation  ( 10)  one  can  obtain  T 
in  term  of  X. 


T  = 


3 


(\-X)(l-aX-bX2  -cX' )' 


(29b) 


Now,  let  us  define 


T'crit  ~  ^X=0.5  ~ 


(1  -  0.5a  -  0.25b  -  0. 125c)" 


(29c) 


Then,  if  T*  <  T*ril,  Equation  (14)  is  used  to  determine  the 
pressure  loss,  otherwise  Equation  (16)  is  employed. 

For  this  MR  fluid  damper’s  geometry,  the  volume  flow 
rate  can  be  expressed  as: 


Q  =  NQm=-V(D1p-D1s)  (30) 

where  N  is  the  number  of  the  MR  valves,  and  Q,„  is  the  vol¬ 
ume  flow  rate  through  each  MR  valve.  Using  Equations  (14) 
and  ( 1 6),  the  pressure  drop  across  one  MR  valve  can  be  writ¬ 
ten  as: 


4xvL, 

"-•“sr* 


In  + 1  8<2, 


V' 


n  jtD, 


A P,„r=- 


i3 

V  "  ) 

4xvL, 


4AL| 

A 


for  T  <Tcril  (31) 
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i-A) 


3«+i  se, w  r  i  r 

n  kD2  j  [Vyj 


for  T  >TC, 


(32) 


The  force  output  is  a  result  of  a  pressure  drop  across  the 
piston  and  can  be  written  as 


F  =  AP-(Dl-D;)  +  F{  sgn(V)  (33) 

4 

where  Ft  is  the  damper  seal  friction  force  (nearly  50  N). 

A  comparison  of  experimental  and  analytical  results  is 
presented  in  Figure  7,  for  a  sinusoidal  motion  with  a  fre¬ 
quency  of  0.5  Hz  and  amplitude  of  1  cm.  The  analytical  re¬ 
sults  are  obtained  directly  from  Equations  (27),  (28),  (3 1)  or 
(32)  and  (33).  Table  3  lists  the  material  and  geometric  prop- 
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Proposed  Simplified 
Hershel-Bulkley  Model 

Gavin  et  al.  (1996) 


ooi  aio  too  looo  100.00  1000.00 

Dimensionless  Yield  Stress,  T* 

Figure  4.  Comparison  between  the  proposed  model  and  Gavin's 
(1996)  approximated  solution  of  Philips  (1969)  Bingham  Plastic 
model,  n  =  1. 


dP.  =  B—+{— — — 1  — ;  for  —<0.5  (22) 
dz  h  v  n  wh2 )  h  h 


h\\-B0 


2n  +  l  2Q  >[  k  f 


n  wh2 )  1% 


71'  t>0'5 


where, 

2«  +  l  3n(\+2n)(\-n)  3nQr+0  ^  B 
«  + 1  16(«  +  1)2  40(n  + 1) 

The  simplified  model  presented  in  Equations  (22)  and  (23) 
provides  a  general  expression  for  the  analysis  of  flow  be¬ 
tween  parallel  plates  that  can  be  employed  as  a  useful  tool  in 
the  design  of  ER  and  MR  devices.  The  Herschel-Bulkley 
model  developed  in  this  study  can  be  reduced  to  the  Bingham 
plastic  model  by  setting  n  =  1.  Using  the  values  of  Table  2,  p 
=  2.629  and  B'  =  3.952,  and.  Equation  (22)  and  (23),  can  be 
respectively  reduced  to: 

dp  _\2kQ+2^5h.  for  0.5  r>A\ 


dz  wh 3 


h\  1-0.593 


K 

for  -2- >0.5  (25) 
h 


the  exact  solution  of  the  cubic  Bingham  flow  equation,  with  a 
relative  error  of  ±2%  for  the  range  0<T*  <1000.  Here,  Equa¬ 
tions  (24)  and  (25)  are  compared  to  those  of  Gavin  et  al. 

( 1996).  For  comparison  purposes,  the  results  from  Equations 
(24)  and  (25)  are  converted  to  the  same  dimensionless  vari¬ 
ables  P  =  (whWUQk)  and  T*  =  (wh2Ty)/(l2Qk)  adopted 
by  Gavin  et  al.  (1996).  Figure  4  shows  excellent  agreement 
between  the  proposed  simplified  Herschel-Bulkley  model 
presented  in  this  study  and  those  of  Gavin  et  al.  (1996).  It 
should  be  noted  that  the  working  range  of  proposed  Equa¬ 
tions  (24)  and  (25)  is  not  limited  to  0  <  T*  <  1000. 

EXPERIMENTAL  VERIFICATION  OF  THE 
HERSCHEL-BULKLEY  MODEL 

The  theoretical  Herschel-Bulkley  model  presented  in  this 
study  is  verified  by  a  comparison  between  the  theoretical  and 
experimental  results.  The  experimental  results  are  obtained 
using  a  new  MR  fluid  damper  developed  at  UNR 
,  (Gordaninejad  and  Breese,  1998).  A  cross-sectional  view  of 
the  UNR  MR  fluid  damper’s  piston  is  shown  in  Figure  5. 
This  figure  shows  the  fluid  flow  path  through  the  piston  as 
well  as  the  geometric  dimensions. 

The  rheological  property  of  MR  fluid  (LORD  132LD)  is 
measured  by  a  universal  dynamic  spectrometer  (Li  et  al., 
1999).  Figure  6  shows  the  shear  stress  versus  shear  strain  rate 
at  various  magnetic  field  strengths.  It  is  found  that  at  a  range 


piiton  Velocity,  V 


Gavin  et  al.  (1996)  presented  an  approximate  formula  to  Figure  5.  Schematic  of  flow  path  for  a  UNR  MR  fluid  damper. 
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Experimental  Results 


Figure  7.  Comparison  between  theoretical  and  experimental 
results  for  a  sinusoidal  motion  at  0.5  Hz  frequency  and  1.0  cm  (0.4 
in.)  amplitude  for  different  input  electric  current  0.0,  1.0,  1.5  and  2.0 
amps. 

erties  used  in  the  theoretical  analysis.  As  demonstrated  in 
Figure  7,  the  Herschel-Bulkley  fluid  model  can  accurately 
predict  the  damping  force  of  the  MR  damper  for  different 
levels  of  input  current. 

SUMMARY  AND  CONCLUSIONS 

The  rheological  behavior  of  ER  and  MR  fluids  can  be 
identified  as  two  distinct  pre-  and  post-yield  regimes.  In  most 
ER  and  MR  devices  such  as  dampers  the  dominant  fluid  be¬ 
havior  occurs  in  the  post-yield  regime.  ER  and  MR  materials 
under  shear  thinning  or  shear  thickening  effects  exhibit  a 
non-linear  behavior  in  the  post-yield  regime.  To  predict  the 
flow  rate  of  these  fluids,  an  accurate  constitutive  model  is 
needed  in  which  the  behavior  of  the  fluid  under  shear  thin¬ 
ning  or  shear  thickening  is  taken  into  consideration. 

In  this  study,  a  detailed  analysis  of  ER  and  MR  field-con¬ 
trollable  fluid  flow  through  pipes  and  parallel  plates  was  pre¬ 
sented  using  the  Herschel-Bulkley  model.  A  steady,  one-di- 
mensional  laminar  flow  was  assumed.  Dimensionless 
variables  were  used  to  simplify  the  continuity  and  momen- 


Table  3.  Material  and  geometric  properties  used  in  the 
Herschel-Bulkley  theoretical  analysis. 


k 

430  Pa  sn 

n 

0.4 

N 

4 

M 

0.4  Pas 

Dp 

0.0508  m 

Ds 

0.0127  m 

D, 

0.00635  m 

d2 

0.00635  m 

J-1 

0.00292  m 

1-2 

0.0939  m 

turn  equations.  The  non-dimensional  plug  thickness  was  de¬ 
termined  to  be  a  crucial  parameter  in  developing  the  simpli¬ 
fied  expressions  for  the  analytical  solution.  When  the  non-di¬ 
mensional  plug  thickness  is  less  than  0.5,  linear  asymptotical 
forms  of  non-dimensional  equations  were  developed.  When 
the  non-dimensional  plug  thickness  is  more  than  0.5,  smaller 
volumetric  flow  rates  can  be  expected.  In  this  case,  a  simpli¬ 
fied  exponential  form  was  derived. 

It  was  demonstrated  that  the  Herschel-Bulkley  model  in 
this  study  can  be  reduced  to  Bingham  plastic,  power  law,  and 
Newtonian  fluid  models,  by  setting  yield  stress,  Tv,  and  in¬ 
dex,  /i,  to  certain  specific-case  values.  In  addition,  the  pro¬ 
posed  simplified  Herschel-Bulkley  model  for  steady  flow  in 
pipes  and  parallel  plates  can  provide  a  convenient  and  gener¬ 
alized  mathematical  form  for  analysis  and  design  of  ER  and 
MR  devices. 
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ABSTRACT 

This  work  presents  a  theoretical  model  for  the  damping  force  of  a  magneto-rheological  fluid  (MRF)  shock 
absorber  of  an  off-road  motorcycle.  The  Bingham  plastic  model  and  a  three-dimensional  electromagnetic 
finite-element  analysis  are  employed  to  develop  a  theoretical  model  to  estimate  the  damping  force  of  a 
MRF  shock  absorber.  The  model  is  based  on  the  physical  parameters  of  the  device  as  well  as  the 
properties  of  the  fluid,  making  it  a  valuable  tool  in  shock  absorber  design  for  a  particular  application.  By 
comparing  the  theoretical  and  experimental  results,  it  is  demonstrated  that  the  model  accurately  predicts 
the  damping  force. 

Keywords:  Magneto-rheological  fluids,  semi-active,  field-controllable,  damper,  off-road  motorcycle,  suspension 


1.  INTRODUCTION 

As  a  component  of  a  semi-active  suspension  system,  the  role  of  the  controllable  shock  absorber  is  to  have  the 
capacity  to  produce  variable  damping  proportionally  reactive  to  the  input  motion,  which  is  often  a  random 
excitation.  Semi-active,  controllable  dampers  are  being  increasingly  investigated  for  various  applications1, .  One 
class  of  controllable,  semi-active  dampers  utilizes  a  magneto-rheological  fluid  (MRF)  as  the  working  fluid.  MRFs 
are  controllable  fluids  that  typically  consist  of  a  silicone-based  carrier  fluid  with  small  iron-based  particles  in 
random  suspension3.  The  application  of  a  magnetic  field  to  the  MRF  causes  the  particles  to  align  themselves  as 
r-hains  along  the  lines  of  magnetic  flux.  The  formation  of  these  chains  causes  a  resistance  to  fluid  flow  due  to  a 
higher  fluid  shear  stress  developed  at  the  wall.  The  higher  fluid  shear  stress  due  to  the  application  of  the  magnetic 
field  is  referred  to  as  a  “MRF  valve”4,5.  The  MRF  valve  can  be  controlled  by  varying  the  input  electric  current  to  the 
device’s  built-in  electromagnet 
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The  objective  of  this  research  is  to  develop  a  theoretical  model  for  a  proof-of-concept  effort.  A  fluid  mechanics 
model  has  been  developed  based  on  the  materials  and  geometric  properties  of  the  device  in  order  to  predict  damping 
characteristics.  The  experimental  component  of  this  study  aids  in  the  validation  and  effectiveness  of  the  developed 
theoretical  model.  This  paper  focuses  on  the  theoretical  aspect  of  this  research. 

The  theoretical  model  is  based  on  the  Bingham  plastic  theoiy  to  develop  the  pressure  drop  across  the  piston  for  the 
case  where  the  MRF  is  in  flow  mode.  In  order  to  determine  the  effect  of  the  magnetic  field  on  the  MRF,  a  three- 
dimensional  finite  element  analysis  (FEA)  is  employed.  The  FEA  results  provide  a  function  that  describes  the 
magnetic  saturation  of  the  fluid  as  a  function  of  input  electric  current.  Integrating  the  electromagnetic  analysis  with 
the  fluid  mechanics  analysis  yields  a  theoretical  model  based  on  the  physical  parameters  of  the  individual  device,  the 
properties  of  the  fluid,  and  the  input  electric  current.  Furthermore,  the  fluid  mechanics  model  accounts  for  the 
external  accumulator  and  piston  by-pass  valving  on  the  compression  stroke.  A  full-scale  MRF  shock  absorber  for 
the  rear  suspension  of  a  motorcycle  is  designed,  manufactured,  and  tested  based  on  a  University  of  Nevada,  Reno’s 
MRF  damper  design6’7.  This  damper  serves  as  the  testbed  for  the  development  of  the  fluid  mechanics  model. 


2.  THEORETICAL  MODELING 

The  pressure  drop  across  the  piston  can  be  determined  using  Bernoulli's  equation  with  an  added  term  to  account  for 
the  fluid  acceleration8,9: 

1  9  It  ldV 

+  &\  =P2  +  g%2+  pHeadLosses  +  p\ — dX  (1) 

2  2  q  at 

Referring  to  Figure  1,  states  1  and  2  have  the  same  fluid  velocity,  and  the  change  in  gravity  between  the  two  points 
can  be  neglected,  so  the  pressure  drop  across  the  piston  can  be  simplified  as: 

LdV 

A P = p(HeadLosse8r  J — dX)  (2) 

o  dt 

Since  either  side  of  the  piston  has  different  areas,  the  damping  force  is  determined  by  taking  the  difference  between 
the  products  of  the  pressure  and  area  on  the  front  and  back  sides  of  the  piston. 

The  flow  channels  in  the  MRF  motorcycle  damper  are  relatively  complex.  Because  it  is  sometimes  difficult  to 
determine  the  overall  pressure  drop  developed  across  the  piston,  it  may  be  helpful  to  examine  the  electric  circuit 
analogue,  where  A P  is  voltage,  volumetric  flow  rate  is  the  electric  current,  and  the  head  losses  are  resistors.  The 
procedure  is  as  follows: 

1 .  Determine  each  head  loss. 

2.  Factor  out  a  local  volumetric  flow  rate  term  [QiocaiKViocai)(Area)=Current],  where  the  coefficients  of 
each  Qiocai  are  the  resistors. 

3.  Combine  the  parallel  and  series  resistors  to  form  a  R^  (equivalent  circuit  with  only  one  resistor  and 
one  volumetric  flow  rate). 

4.  Determine  the  pressure  drop,  AP=(p)(Req)(Qlotol),  where  QtouJ  -  overall  conserved  volumetric  flow  rate. 
A  graphical  representation  of  this  scheme  is  presented  in  Figure  2. 


The  theoretical  model  presented  is  developed  to  predict  the  behavior  of  a  new  prototype  MRF  damper  (MRF-MC1) 
for  the  rear  suspension  of  a  motorcycle.  A  cross-sectional  view  of  the  MRF  motorcycle  damper  is  shown  in  Figure 
1.  The  fluid  flow  path  through  the  piston  as  well  as  the  critical  parametric  variables  used  in  the  theoretical 
formulation  are  demonstrated  in  the  Figure.  Referring  to  Figure  1,  the  MRF-MC1  damper  has  an  external 
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accumulator  with  a  pressure-filled  bag  that  defonns  to  absorb  the  added  volume  to  the  system  due  to  the  rod.  Also, 
the  damper  has  piston  by-pass  valving  on  the  compression  stroke.  These  two  features  must  be  addressed  in  the 
development  of  die  model.  The  first  step  in  the  derivation  of  the  model  is  to  find  the  pressure  drop  across  the  piston. 
Then,  die  effects  of  the  accumulator  are  integrated  into  the  model  to  find  the  overall  force  of  the  damper. 

Due  to  the  effects  of  by-pass  valving  on  the  compression  stroke,  the  pressure  drops  across  the  piston  are  different  on 
the  compression  and  rebound  strokes.  As  the  wave  washer  on  the  piston  opens  (corresponding  to  large  pressure 
gradients),  the  diameter  of  each  by-pass  valve  changes  from  0  to  D3.  This  change  is  proportional  to  the  velocity 
which  means  that  when  the  velocity  is  very  high,  the  wave  washer  may  be  fully  open.  A  function  is  selected  to 
describe  the  behavior  of  the  wave  washers  under  a  sinousoidal  input  motion. 


P  ~  P  0  cos(fitf)  (3) 

po  describes  how  open  each  by-pass  valve  is,  P  ranges  from  0  to  po.  When  P=0,  the  valve  is  fully  closed.  This 
1  would  correspond  to  the  rebound  stroke.  When  p=po,  the  valve  opening  is  a  maximum  at  p0  (ie.:  compression 

stroke).  The  by-pass  valve  diameter  is  equal  to  PD3. 

There  are  six  resistances  to  flow  through  the  piston  of  the  MRF-MC1  damper  on  the  compression  stroke8,9.  These 
resistances  are  due  to  frictional  effects,  inertial  effects,  and  fluid  bifurcation  effects.  The  piston  resistance  network 
is  shown  graphically  in  Figure  3.  Combining  the  series  and  parallel  a  single  equivalent  resistance  (Rcq)  is  obtained. 
I  Multiplying  the  equivalent  resistance  by  the  overall  conserved  volumetric  flow  rate  yields  the  pressure  drop  across 

the  piston  as: 


32M,(P,‘-.P,1h  32  N^DX 

0®,)4  C(J3D,y 

where  Vt  is  the  local  fluid  velocity  through  the  center  of  the  piston,  and  given  as10: 


VfjJK  -g/h  32^C0(C |J2  -Kygggv 


(WDj)* _ £> _ 

64I|fe  |  C14  MoPt2  ,  32A( |  32Afr„4I>|2 

D,2  D  D>  D*  O/Bj)4 


(4) 


(5) 


Some  of  the  parametric  variables  shown  in  Equations  (4)  and  (5)  are  defined  in  Figure  1.  In  these  Equations,  m  is 
the  zero-field  fluid  viscosity,  £  is  the  number  of  by-pass  valve  ports,  and  N  is  the  number  of  MRF  channel  flow 
ports.  Referring  to  the  second  term  in  the  numerator  of  Equation  (5),  C0  is  a  material  constant.  It  represents  the 
slope  of  the  line  describing  the  shear  stress  in  the  fluid  as  a  function  of  magnetic  field  saturation.  The  term 
(CiI2+C2I)  represents  the  magentic  behavior  of  the  MRF  due  to  input  electric  current. 

The  pressure  backed  accumulator  bag  acts  similar  to  a  spring,  where  the  force  on  the  bag  is  proportional  to  its 
deformation.  Therefore,  it  is  necessary  to  relate  the  deformation  of  the  accumulator  bag  to  the  known  main  piston 
displacement.  Volumetric  flow  conservation  allows  this  relation  to  be  developed. 

Because  of  the  accumulator,  there  are  two  different  volumetric  flow  rates  that  are  conserved  in  the  damper.  The  first 
is  the  flow  that  is  conserved  through  the  piston.  The  second  is  a  flow  generated  by  the  increased  volume  to  the 
system  due  to  the  piston  rod,  and  is  given  as: 
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The  flow  due  to  the  increase  in  shaft  volume  is  sent  to  the  accumulator.  Equation  (6)  shows  that  die  amount  of  fluid 
sent  to  the  accumulator  is  a  function  of  the  main  piston  displacement. 


The  accumulator  can  be  approximated  as  an  isentropic  process,  meaning  that  due  to  minimal  friction,  the  system  is 
internally  reversible  and  that  the  system  is  nearly  adiabatic".  Equation  (7)  presents  the  pressure  in  the  accumulator 
as  a  function  of  main  piston  displacement  (Xp): 


Pa  (0  —  Pch  arg 


(7) 


where  Pcha^e  and  VChargt  correspond  to  the  initial  pressure  and  and  volume  of  the  accumulator  bladder,  respectively. 
It  should  be  noted  that  the  piston  is  fully  extended,  and  X  corresponds  to  the  displacement  of  the  main  piston  from 
its  fully  extended  position  to  its  oscilation  midpoint. 


To  determine  the  force  on  the  piston,  the  actual  pressures  on  the  front  and  back  of  the  piston  must  be  multiplied  by 
the  corresponding  cross-sectional  areas.  This  gives  a  force  on  the  front  and  a  force  on  the  back  of  the  piston.  The 
actual  force  on  the  piston  is  the  difference  of  the  two.  In  other  words: 


where: 


f piston  P\^p\  PlApl 


4»  =fk2-A2] 


a  n  2 

Ap2~4Dp 


Equation  (9)  can  be  expanded  as: 


F piston  piston 


Ij-k+AP.^, 


In  expanded  form,  the  force  on  the  piston  is: 


(8) 

(9) 

(10) 


(11) 
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where  K9  and  K,  |  represent  dimensionless  entrance  and  discharge  loss  cnpffirWc  a 

cylinder  to  the  accumulator,  respectively.  8  °S  coefflcients  of  *e  Port  connectmg  the  main 


3.  RESULTS 

ssi-*  ■  *”  “rycfe  *■*-  <MRr-Mc» 

show,  4.  d*n*««d  «  ,h.2k  **?«■  ™<  5 

many  tests  to  explore  the  effects  of  input  stroke  length  inmit  ^  exPenmental  damper  underwent 

the  force  generated  by  die  unit.  The  results  of  the  teste  are  for  ^  a-  T’  ^  mput  cuirent  *°  electromagnet  on 
done  by  the  damper.  Figures  6  and  7  show  the  comnari™t  secernent  curves,  which  represents  the  work 

displacement  loops  as  a  function  of  input  current  at  an  input  stroke  of  experimental  force' 

Hz,  respectively.  Z  ?£T“ 1J 

(0.00254  m>,  L,  -  0.65  in.  ,0.0165  "S’ U  'V  * 

(861.8  kPa),  p  =  3004  kg/(mA3),  WD  =  60.0,  K3=05L=10K  -  inr  i  n  J ’  ?E  5  m  ’  Pcharge  125  psi 
Oersted/ amp2,  C  2  =  97.7  Oersted/^p,  N  =  4,  £  6  *  ’  **  ~  L°’ K"  =  10>  C*  =  *  Pa/Oersted,  C,  =  32.7 


Both  Figures  6  and  7  show  excellent  agreement  between  the  the™-„t,v„i  j 

differences  between  the  theoretical  and  experimental  results  at  the  ■  d  expefmental  results-  The  minor 
amounts  of  air  in  the  damper.  expenmental  results  at  the  maximum  stroke  lengths  are  due  to  small 


for  fc  re„  ^ion  of 

analysis  was  performed  to  incorporate  the  magnetic  effects  nf  mt?f  i  •°iy'  An  eiectromagnetic  finite  element 
force  is  estimated  for  the  dime*  ions  of  the°retiCal  modeL  ’^.‘lamping 

proposed  model  is  checked  by  a  comparison  between  the  theoretical  and  e  103  m0t,°r  CycIe'  7516  validity  of  the 
exceUent  agreement.  theoretical  and  experimental  results,  which  demonstrated 
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Figure  1.  Schematic  of  flow  path  for  MRF-MC1  damper  with  by-pass  valving. 


(Qlocal)2  ^  (Qlocal)2 


Figure  2.  Graphical  representation  of  the  electric  circuit  analogue  procedure. 


Figure  3.  Graphical  representation  of  the  compression  stroke  resistance  network 


Figure  4.  Basic  design  and  dimensions  of  MRF-MC1  damper. 


Figure  5.  MRF-MC1  shock  absorber  mounted  on  Instron  shock  dynamometer. 


Experimental 


Figure  6.  Comparisons  between  theoretical  and  experimental  results 
inusoidal  motion  with  1.5  Hz  frequency,  and  0.508  cm  amplitude  mput  for 
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Abstract 


This  experimental  study  presents  the  development 
md  evaluation  of  a  conlrollable,  semi-active  magneto- 
rheological  fluid  (MRF)  shock  absorber  for  a  rear 
suspension  of  an  off-road  motorcycle.  The  MRF  device 
is  designed  and  constructed  at  the  University  of  Nevada, 
Reno  (UNR).  An  experimental  study  ts :  conducted  to 
examine  the  performance  of  the  proposed  retrofit  MRF 
damper.  The  experimental  results  show  that  the  MRF 
damper  can  provide  wide  range  damping  forces. 


INTRODUCTION 

As  a  component  of  a  semi-active  suspension  system, 
the  role  of  the  controllable  shock  absorber  is  to  have  the 
capacity  to  produce  variable  damping  force 
proportionally  reactive  to  an  input  excitation.  Semi¬ 
active,  controllable  dampers  are  bemg  investigated 
increasingly  for  various  applications.  One  class  of 
controllable,  semi-active  dampers  is  MRF  shock: 
absorbers.  Magneto-rheological  fluids  typically  consist 
of  a  silicone-based  carrier  fluid  with  small  iron-based 
particles  in  random  suspension  [1],  The  application  of  a 
magnetic  field  to  the  MRF  causes  the  particles  to  ahgn 
as  chains  along  the  lines  of  magnetic  flux.  The 
formation  of  these  chains  causes  a  resistance  to  fluid 
flow  due  to  the  development  of  higher  fluid  shear 
stresses.  By  changing  the  magnetic  field,  the  shear 
yield  stress  of  the  fluid  changes;  thus  producing  a 
controllable  force  for  the  MRF  damper.  In  other  words, 
it  is  possible  to  control  the  damping  forces  of  a  device 
utilizing  MRF,  by  varying  the  input  electric  current  to 
the  device’s  built-in  electromagnet. 

The  ability  to  vary  damping  for  a  wide  range  of 
inputs  is  the  main  feature  of  semi-active  suspension 
systems.  Other  features  include  control  of  vehicle  pitch 
during  braking  and  acceleration  and  accommodation  for 
a  variable  vehicle  mass  (i.e.,  light  or  heavy  riders).  In 
comparison  to  typical,  passive  suspension  systems, 
semi-active  shock  absorbers  allow  for  softer  damping 
when  needed,  in  addition  to  harder  damping  for 
situations  that  demand  it  The  main  benefits  of  the  MR 
technology  for  off-road  motorcycles  are:  improved 
operator  ride  comfort,  an  increased  level  of  vehicle 


control  due  to  a  lower  disturbance  input,  and  faster 
operational  speeds  over  similar  terrain. 

Semi-active  vibration  control  devices  are  receiving 
significant  attention  because  they  can  offer  combined 
advantages  of  both  passive  and  active  systems. 
Examples  of  such  devices  include  variable  orifice 
dampers,  variable  stiffness  devices,  and  friction 
controllable  isolators  and  braces.  A  class  of  semi-active 
controllable  devices  utilizing  MRF  has  been  established. 
There  have  been  numerous  studies  on  MR  fluids  and 
devices  in  recent  years  due  to  the  realization  that  there 
are  many  potential  applications  [1-8]. 

A  new  generation  of  MRF  dampers  has  been 
designed  and  developed  at  UNR  that  provide  robust, 
cost-effective,  vibration  damping  solutions.  Several 
innovative  prototype  MRF  dampers  have  been 
developed  and  examined  at  the  UNR’s  Composite  and 
Intelligent  Materials  Laboratory  for  different 
applications,  such  as  mountain  bikes,  motorcycles,  off¬ 
road  vehicles,  and  bridges  [5-8]. 

This  study  focuses  on  a  proof-of-concept 
experimental  study  to  demonstrate  the  capabilities  of 
UNR  MRF  retrofit  damper  for  an  off-road  motorcycle. 
A  full-scale  MRF  retrofit  shock  absorber  for  the  rear 
suspension  of  a  motorcycle  is  designed,  manufactured, 
and  examined. 

EXPERIMENTAL  STUDY 

ai  DESIGN  METHODOLOGYl  In  this  study,  a 
specific  design  methodology  is  followed  in  the 
development  of  the  MRF  retrofit  damper.  This  design 
methodology  is  based  on  the  fail-safe  characteristic  of 
the  MRF  devices.  MRF  dampers  are  fail-safe  in  their 
operation,  meaning  that  in  zero-field  (off-state),  passive 
mode  operation  (i.e.,  no  electric  current  supply  to  the 
damper)  the  MRF  shock  absorber  performs  similar  to 
the  original  equipment  manufacture  (OEM)  shock 
absorber.  Therefore,  the  steps  in  designing  a  MRF 
retrofit  shock  absorber  for  the  rear  suspension  of  an  off¬ 
road  motorcycle  are: 

1.  Benchmark  test  the  existing  OEM  shock  absorber. 

2.  Design  MRF  damper  to  give  similar  performance  at 
zero-field  as  the  OEM  unit. 


267 


3.  Apply  power  to  the  MRF  damper  to  generate  an 

“envelope”  of  increased  performance,  with  the 

OEM  performance  near  the  bottom  of  the  envelope. 

An  OEM  rear  suspension  shock  absorber  of  an  off¬ 
road  motorcycle  is  characterized  in  order  to  determine 
the  base  performance  of  the  unit  as  a  function  of  input 
amplitude  of  the  motion  and  frequency.  After  fully 
characterizing  the  OEM  damper,  it  is  disassembled  and 
carefully  studied.  The  OEM  piston  unit  is  removed 
from  the  OEM  shock  absorber  body  and  a  MRF  retrofit 
piston  is  assembled  in  its  place.  Replacing  the 
hydraulic  fluid  with  MR  fluid  completes  the  retrofit  of 
the  OEM  damper.  The  result  is  a  shock  absorber  that 
has  the  same  dimensions  as  the  OEM  damper  with  a 
MRF  piston  retrofit  Figure  1  shows  the  components  of 
the  OEM  shock  absorber,  and  Figure  2  shows  the 
components  of  the  MRF  retrofit  damper. 

b)  MRF  SHOCK  ABSORBER  DESIGN:  The  design 
of  the  MRF  shock  absorber  involves  utilizing  a 
parametric  non-Newtonian  fluid  mechanics  model  and 
extensive  electro-magnetic  modeling  by  means  of  a 
three-dimensional  magnetic  finite  element  analysis 
(FEA).  For  simplicity  in  manufacturing,  and  also  to 
demonstrate  the  feasibility  of  retrofitting  an  OEM  shock 
absorber  with  a  MRF  piston,  the  OEM  shock  absorber, 
accumulator  bag,  seal  pack,  and  chrome-plated  rod  are 
used  in  the  construction  of  the  UNR  MRF  retrofit 
damper.  It  is  important  to  preserve  the  3 -inch  stroke  of 
the  damper,  which  means  that  the  MRF  piston  must  be, 
at  most,  the  same  size  as  the  OEM  piston.  The 
following  are  the  physical  design  constraints  for  the 
MRF  piston: 

1 .  Maximum  length  =  4.064  cm  ( 1 .6  inches) 

2.  Maximum  diameter  =  4.3 1 8  cm  ( 1 .7  inches) 

3.  Piston  ring  grove  width  =  0.8128  cm  (0.32  inches) 

Furthermore,  the  MRF  piston  retrofit  must  have 
compression  by-pass  valving  to  protect  the  damper  from 
harmful  impulse  loads  on  the  compression  stroke.  With 
the  necessary  physical  constraints  of  the  piston  clearly 
defined  as  shown  above,  the  next  step  in  the  design  is  to 
determine  the  dimension  of  the  fluid  flow  ports  to  give 
damping  forces  that  are  similar  to  the  OEM  shock 
absorber.  A  theoretical  non-Newtonian  fluid  mechanics 
model,  developed  by  the  authors,  is  employed  to 
estimate  the  size  of  the  flow  channels  [7].  The  model 
predicts  the  damping  force  based  on  the  physical 
parameters  of  the  fluid  and  the  device.  Using  this 
model,  the  dimensions  of  the  fluid  flow  channels  in  the 
MRF  piston  can  be  determined  so  that  the  MRF  retrofit 
shock  absorber  produces  similar  damping  forces  as  the 
OEM  damper  at  zero-field  state. 

With  the  physical  parameters  of  the  MRF  shock 
absorber  specified,  the  final  step  in  the  design  is  to  use 
FEA  to  design  and  characterize  the  piston’s  built-in 


electromagnet.  In  particular,  the  effects  of  material 
properties  on  the  magnetic  flux  lines  and  magnetic 
saturation  of  the  MRF  is  investigated  by  FEA.  In 
addition,  by  varying  the  electric  input  current  in  the 
FEA  model,  a  function  is  generated  describing  the 
shear  yiels  stress  of  the  MRF  as  a  function  of  input 
electric  current.  This  provides  a  controllable,  current- 
dependent  term  in  the  fluid  mechanics  model  which 
accurately  predicts  the  current  dependent  variability  of 
the  MRF  damper.  Upon  completion  of  the  fininte 
element  analysis,  the  most  appropriate  material 
configuration  for  the  construction  of  the  MRF  damper, 
is  determined. 

By  both,  non-Newtonian  fluid  mechanics  modeling 
and  electro-magnetic  finite  element  analysis,  a  MRF 
retrofit  piston  is  designed.  The  design  preserves  the 
OEM  stroke  length  and  performance  and  incorporates 
compression  stroke  by-pass  valving  to  protect  the 
device  from  harmful  impulse  loading.  Next,  a  prototype 
MRF  piston  is  manufactured  and  characterized. 

r)  MRF  DAMPER  CHARA  CTERIZA  TIQM  Figure  3 
shows  a  photo  of  the  experimental  set  up  with  UNR 
MRF  retrofit  shock  absorber  mounted  on  the  table.  The 
shock  absorbers  for  this  study  are  tested  on  an  Instron 
8821  S  hydraulic  shock  testing  dynamometer.  The 
Instron  has  a  maximum  stroke  of  15.24  cm  (6.0  inches), 
and  it  is  equipped  with  a  22  KN  (5,000  Ibf)  load  cell 
that  measures  the  damping  forces  of  the  shock  absorber 
and  a  displacement/velocity  transducer  measures  the 
displacement  and  velocity  of  the  actuator.  Experimental 
data  is  collected  at  a  sampling  rate  of  200  Hz. 

Using  the  Instron  shock  dynamometer,  each  of  the 
shock  absorbers  evaluated  in  this  work  is  undergone 
various  performance  tests.  A  performance  test  consists 
of  a  frequency  sweep  held  at  constant  displacement 
input  Each  specific  amplitude  and  frequency  represents 
a  run  of  the  performance  test  The  results  of  a 
performance  test  are  a  force-displacement  and  force- 
velocity  response  curves.  These  curves  represent  die 
work  done  and  the  energy  absorbed  by  the  damper, 
respectively. 

RESULTS  AND  DISCUSSION 

The  performance  of  the  UNR  MRF  retrofit  damper 
for  rear  suspension  of  an  off-road  motorcycle  is 
evaluated  by  conducting  numerous  tests.  The 
experiments  are  mainly  performed  to  assess  the  effects 
of  input  amplitude,  input  frequency,  and  input  electric 
current  to  the  electromagnet,  on  the  damping  force 
generated  by  the  MRF  damper.  The  results  of  these 
tests  are  a  force-displacement  performance  curve  that 
represents  the  work  done  by  the  damper. 

In  each  of  the  force-displacement  loops  presented  © 
this  work,  the  positive  force  and  displacement  values 
correspond  to  the  rebound  stroke  whereas  the  negaO« 
force  and  displacement  values  correspond  to 
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Lcompression  stroke.  The  compression  stroke  force 
values  do  not  increase  in  value  with  increasing 
frequency  as  much  as  the  rebound  force  values.  This  is 
one  effect  of  by-pass  valving  on  the  compression  stroke. 

As  the  input  frequency  increases  (corresponding  to  large 
piston  velocities  and  pressure  drops),  by-pass  valves 
open  to  reduce  the  pressure.  The  by-pass  valves  on  the 
rebound  stroke  are  not  as  dominant,  so  force  grows 
much  higher  with  increasing  frequency. 

Figure  4  shows  the  force-displacement  loops  for  the 
MRF  retrofit  shock  absorber  as  a  function  of  input 
electric  current  This  plot  demonstrates  that  the 
maximum  rebound  damping  force  changes  from  about 
30  N  to  430  N,  between  0  and  2  Amp  input  electric 
current.  In  other  words,  the  “performance  envelope  is 
1  represented  by  a  400  N  window  of  controllability. 

Figure  5  shows  the  force-displacement  comparison 
between  the  OEM  damper,  the  UNR  MRF  retrofit 
damper  at  zero-field,  and  the  MRF  damper  with  2  amps 
current  input  at  amplitude  of  2.5  mm  (0.2  in.)  and  1  Hz 
input  frequency.  Note  the  performance  envelope 
created  around  the  OEM  shock  absorber  performance. 

1  Again,  this  envelope  represents  the  controllability  of  the 
MRF  damper.  The  zero-field  (off-state)  force  is  slightly 
lower  than  OEM  force  and  the  force  generated  by  the  2 
Amp  electric  current  input  is  higher  than  the  OEM 

Finally  Figure  6  shows  the  peak  rebound  force 
comparison  between  the  OEM  and  the  MRF  damper  at 
1  zero-field,  and  the  MRF  damper  with  2  Amps  electric 
current  input  as  a  function  of  excitation  frequency.  The 
input  amplitude  for  these  results  is  2.5  mm  (0.2  in.). 
Again,  note  that  the  0  amp  and  2  amp  curves  form  an 
envelope  around  the  OEM  damper  curve.  These  results 
demonstrate  the  design  methodology  discussed 
previously.  The  OEM  damper  is  characterized  as  a 
function  of  input  frequency.  Using  the  fluid  mechanics 
model,  the  MRF  retrofit  damper  is  designed  to  produce 
zero-field  force  values  that  are  slightly  less  than  the 
OEM  damper  independent  of  input  frequency.  With  the 
application  of  2  Amps  input  electric  current  to  me 
damper’s  built-in  electromagnet,  the  MRF  damping 
levels  increase  by  400  N,  creating  a  significant  dynamic 
i  force  range  for  the  MRF  retrofit  damper. 

CONCLUSIONS 

A  UNR  MRF  retrofit  damper  is  designed, 
constructed  and  evaluated  for  an  off-road  motorcycle. 
The  MRF  shock  absorber  has  similar  performance  to 
\  that  of  an  OEM  damper  at  zero-field.  With  the 
application  of  electric  current  to  the  damper  s  built-in 
electromagnet,  a  wide  controllable  range  of  damping 
force  above  the  OEM  performance  is  produced.  The 


feasibility  of  implementing  compact,  MRF  technology 
onto  an  off-road  motorcycle  while  maintaining  the 
OEM  stroke  length  performance,  and  compression  by¬ 
pass  valving  characteristics,  is  successfully 
demonstrated. 
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igure  5.  Force-Displacement  loops  for  the  OEM  and  the  MRF  Damper  at  1.0  Hz  frequency  and  2.5  mm  amplitude 
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Figure  6.  Dynamic  force  range  of  the  MRF  damper  at  2.5  mm  amplitude. 
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